A GAME INTERPRETATION OF THE NEUMANN PROBLEM FOR FULLY 
NONLINEAR PARABOLIC AND ELLIPTIC EQUATIONS 



JEAN-PAUL DANIEL 



Abstract. We provide a deterministic-control-based interpretation for a broad class of fully nonlinear 
parabolic and elliptic PDEs with continuous Neumann boundary conditions in a smooth domain. We 
construct families of two-person games depending on a small parameter e which extend those proposed 
by Kohn and Serfaty I'illl . These new games treat a Neumann boundary condition by introducing some 
specific rules near the boundary. We show that the value function converges, in the viscosity sense, to 
the solution of the PDE as e tends to zero. Moreover, our construction allows us to treat both the 
oblique and the mixed type Dirichlet-Neumann boundary conditions. 

1. Introduction 

In this paper, we propose a deterministic control interpretation, via "two persons repeated games", 
for a broad class of fully nonlinear equations of elliptic or parabolic type with a continuous Neumann 



boundary condition in a smooth (not necessarily bounded) domain. In their seminal paper 21], Kohn and 
Serfaty focused on the one hand on the whole space case in the parabolic setting and on the other hand 
on the Dirichlet problem in the elliptic framework. They construct a monotone and consistent difference 
approximation of the operator from the dynamic programming principle associated to the game. 

Our motivation here is to adapt their approach to the Neumann problem in both settings. Furthermore, 
once this issue is solved, we will see how the oblique or the mixed type Dirichlet-Neumann boundary 
problem can also be treated by this analysis. We consider equations in a domain C M. N having the 
form 

(1.1) -u t + f{t,x,u,Du,D 2 u) = 
or 

(1.2) f{x lU ,Du,D 2 u) + Xu = 0, 

where / is elliptic in the sense that / is monotone in its last variable, subject to the Neumann boundary 
condition 

As in [21], the class of functions / considered is large, including those that are non-monotone in the u 
argument and degenerate in the D 2 u argument. We make the same hypotheses on the continuity, growth, 
and u-dependence of / imposed in [2lj. They are recalled at the end of the section. In the stationary 
setting (|1.2p . we focus on the Neumann problem, solving the equation in a domain f2 with (|1.3[) at dVl. In 
the time-dependent setting (|l.lj) . we address the Cauchy problem, solving the equation with (|1.3j) at dtt 
for t < T and u = g at terminal time t — T. The PDEs and boundary conditions are always interpreted 
in the "viscosity sense" (Section [3] presents a review of this notion) . 

Our games have two opposite players, Helen and Mark, who always make decisions rationally and 
deterministically. The rules depend on the form of the equation, but there is always a small parameter 
e, which governs the spatial step size and (in time-dependent problems) the time step. Helen's goal is to 
optimize her worst-case outcome. When / is independent of u, we shall characterize her value function 



Date: December 14, 2012. 



1 



2 



JEAN-PAUL DANIEL 



u e by the dynamic programming principle. If / depends also on u, the technicality of ours arguments 
requires to introduce a level-set formulation since the uniqueness of the viscosity solution is no longer 
guaranteed. The score U e of Helen now depends on a new parameter z6l. In the parabolic setting, it 
is defined by an induction backward in time given by 

Vz e R, U e (x,z,t) = maxminC/ E (K + Ax,z + Az,t + At), 

pS Ax 

endowed with the final-time condition U £ (x,z,t) = g(x) — z. The max on p, T and the min on Ax 
are given by some constrains depending on the rules of the game and some powers of e. This dynamic 
programming principle is similar to the one given in (2lL Section 2.3]. In that case, our value functions 
u e of interest are defined through the 0-level set of U e with respect to z as the maximal and the minimal 
solutions of U £ (x,z,t) — 0. They satisfy two dynamic programming inequalities (for the details of our 
games and the definition of Helen's value function, see Section [2]). 

Roughly speaking, our main result states that 

limsupi/ is a viscosity subsolution of the PDE, and 

liminf u E is a viscosity supersolution of the PDE. 

For the general theory of viscosity solutions to fully nonlinear equations with Neumann (or oblique) 
boundary condition the reader is referred to (l2l. [3L fl9|. As for the Neumann boundary condition, its 
relaxation in the viscosity sense was first proposed by Lions (2^ . 

Our result is most interesting when the PDE has a comparison principle, i.e. when every subsolution 
must lie below any supersolution. For such equations, we conclude that limu £ exists and is the unique 
viscosity solution of the PDE. In the case when / is continuous in all its variable, there are already 
a lot of comparison and existence results for viscosity solutions of second order parabolic PDEs with 
general Neumann type boundary conditions. We refer for this to 0, [2^, and references therein. 
For homogeneous Neumann conditions, Sato [27] has obtained such a comparison principle for certain 
parabolic PDEs. 

We are interested here in giving a game interpretation for fully nonlinear parabolic and elliptic equa- 
tions with a Neumann condition. Applications of the Neumann condition to deterministic optimal control 
and differential games theory in [22] rely much on a reflection process, the solution of the deterministic 



Skorokhod problem. Its properties in differents situations are studied in many articles such as 28|, IH 13 1 . 
The case of the Neumann problem for the motion by mean curvature was studied by Giga and Liu 17j . 
There, a billiard game was introduced to extend the interpretation made by Kohn and Serfaty (2p| via the 
game of Paul and Carol. It was based on the natural idea that a homogeneous Neumann condition will 
be well-modeled by a reflection on the boundary. Liu also applies this billiard dynamics to study some 
first order Hamilton- Jacobi equations with Neumann or oblique boundary conditions [25]. Nevertheless, 
in our case, if we want to give a billiard interpretation with a bouncing rule which can send the particle 
far from the boundary, we can only manage to solve the homogeneous case. This is not too surprising 
because the reflection across dil is precisely associated to a homogeneous Neumann condition. 

Another approach linked to the Neumann condition is to proceed by penalization on the dynamics. 
For a bounded convex domain, Lions, Menaldi and Sznitman [23] construct a sequence of stochastic 
differential equations with a term in the drift coefficients that strongly penalizes the process from leaving 
the domain. Its solution converges towards a diffusion process which reflects across the boundary with 
respect to the normal vector. Barles and Lions [7] also treat the oblique case by precisely establishing the 
links between some approximated processes and the elliptic operators associated to the original oblique 
stochastic dynamics. 

Instead of a billiard, our approach here proceeds by a suitable penalization on the dynamics depending 
on the Neumann boundary condition. It will be favorable to one player or the other according to its 
sign. We modify the rules of the game only in a small neighborhood of the boundary. The particle driven 
by the players can leave the domain but then it is projected within. This particular move, combined 
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with a proper weight associated to the Neumann boundary condition, gives the required penalization. 
Outside this region, the usual rules are conserved. Therefore the previous analysis within £1 done by 
Kohn and Serfaty can be preserved. We focus all along this article on the changes near the boundary and 
their consequences on the global convergence theorem. In this context, the modification of the rules of 
the original game introduces many additional difficulties intervening at the different steps of the proof. 
Most of all, they are due to the geometry of the domain or the distance to the boundary. As a result, 
our games seem like a natural adaptation of the games proposed by Kohn and Serfaty by permitting to 
solve an inhomogeneous Neumann condition h depending on x on the boundary. We only require h to be 
continuous and uniformly bounded, the domain to be C 2 and to satisfy some natural geometric conditions 
in order to ensure the well-posedness of our games. Moreover our approach can easily be extended both 
to the oblique and the mixed Neumann-Dirichlet boundary conditions in both parabolic and elliptic 
settings. Our games can be compared to those proposed in [21| for the elliptic Dirichlet problem: if the 
particle crosses the boundary, the game is immediately stopped and Helen receives a bonus b{xp) where 
b corresponds to the Dirichlet boundary condition and xf is the final position. Meanwhile, our games 
cannot stop unexpectedly, no matter the boundary is crossed or not. 

Our games, like the ones proposed by Kohn and Serfaty, are deterministic but closely related to a 
recently developed stochastic representation due to Cheridito, Soner, Touzi and Victoir [ll| (their work 
uses a backward stochastic differential equation, BSDE, whose structure depends on the form of the 
equation). 

Another interpretation is to look our games as a numerical scheme whose solution is an approximation 
of a solution of a certain PDE. This aspect is classical and has already been exploited in several contexts. 
We mention the work of Peres, Schramm, Sheffield and Wilson [26] who showed that the infinity Laplace 
equation describes the continuum limit of the value function of a two-player, random-turn game called 
e-step tug-of-war. In related work, Armstrong, Smart and Sommersille [2| obtained existence, uniqueness 
and stability results for an infinity Laplace equation with mixed Dirichlet-Neumann boundary terms by 
comparing solutions of the PDE to subsolutions and supersolutions of a certain finite difference scheme, 
by following a previous work of Armstrong and Smart for the Dirichlet case 

This paper is organized as follows: 

• Section [5] presents the two-person games that we associate with the PDEs and (|1.2[) , mo- 
tivating and stating our main results. The section starts with a simple case before adressing 
the general one. Understanding our games is still easy, though the technicality of our proofs 
is increased. Since / depends on u, the game determines a pair of value functions u 6 and v E . 
Section 12.2.11 gives a formal argument linking the principle of dynamic programming to the PDE 
in the limit s — >• and giving the optimal strategies for Helen that will be essential to obtain 
consistency at Section HJ 

• Section [3] addresses the link between our game and the PDE with full rigor. The proofs of 
convergence follow the background method of Barles and Souganidis [10), i.e. they use the 
stability, monotonicity and consistency of the schemes provided by our games. Their theorem 
states that if a numerical scheme is monotone, stable, and consistent, then the associated "lower 
semi-relaxed limit" is a viscosity supersolution and the associated "upper semi-relaxed limit" is 
a viscosity subsolution. The main result in Section [3] is a specialization of their theorem in our 
framework: if v e and u e remain bounded as e — > then the lower relaxed semi-limit of if is 
a viscosity supersolution and the upper relaxed semi-limit of u £ is a viscosity subsolution. We 
also have v e < u e with no extra hypothesis in the parabolic setting, or if / is monotone in u 
in the elliptic setting. If the PDE has a comparison principle (see then it follows that 
limu 6 = limi/ exists and is the unique viscosity solution of the PDE. 

• The analysis in Section [3J shows that consistency and stability imply convergence. Sections [4] 
and [5] provide the required consistency and stability results. The new difficulties due to the 
penalization corresponding to the Neumann condition arise here. The main difficulty is to control 
the degeneration of the consistency estimate obtained in (2l| with respect to the penalization. 
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Therefore we will mainly focus on the consistency estimates whereas the needed changes for 
stability will be simply indicated. 
• Section |6] describes the games associated on the one hand to the oblique problem in the para- 
bolic setting and on the other hand to the mixed type Dirichlet-Neumann boundary conditions 
in the elliptic framework. By combining the results associated to the game associated to the 
Neumann problem in Section [2] with the ideas already presented in |2l| . we can obtain the results 
of convergence. 

Notation: The term domain will be reserved for a nonempty, connected, and open subset of M. N . If 
x, y G R N , (x, y) denotes the usual Euclidean inner product and ||a:|| the Euclidean length of x. If A is 
a N x N matrix, \\A\\ denotes the operator norm \\A\\ = sup ||Ax||. § N denotes the set of symmetric 



N x N matrices and Eij the (i,j')-th matrix unit, the matrix whose only nonzero element is equal to 1 
and occupies the (i,j)-th position. 

Let be a domain in R w and (0) be the vector space of fc-times continuously differentiable functions 
u: — > R, such that all the partial derivatives of u up to order k are bounded on 0. For a domain f2, we 
define 



If £1 is a smooth domain, say C 2 , the distance function to dVl is denoted by d = d(-, dil), and we recall 
that, for all x G <9f2, the outward normal n(x) to dfl at x is given by n(x) = ~Dd{x). 

Observe that, if 80, is assumed to be bounded and at least of class C 2 , any x G M. N lying in a sufficiently 
small neighborhood of the boundary admits a unique projection onto <9f2, denoted by 



In particular, the vector x — x is parallel to n(x). The projection onto fl will be denoted by proj^-. When 
it is well-defined, it can be decomposed as 



For each a > 0, we define il(a) — {x g fl, d(x) < a}. We recall the following classical geometric condition 
(see e.g. [3). 

Definition 1.1 (Interior ball condition). The domain f2 satisfies the interior ball condition at xq 6 dQ 
if there exists an open ball B <Z VL with Xq G dB. 

We close this introduction by listing our main hypotheses on the form of the PDE. First of all we 
precise some hypotheses on the domain fl. Throughout this article, fl will denote a C 2 -domain. In the 
unbounded case, we impose the following slightly stronger condition than the interior ball condition. 

Definition 1.2 (Uniform interior/exterior ball condition). The domain Q satisfies the uniform interior 
ball condition if there exists r > such that for all x € dQ there exists an open ball JJcSl with x £ dB 
and radius r. Moreover, the domain Q satisfies the uniform exterior ball condition ifM. N \tt satisfies the 
uniform interior ball condition. 

We observe that the uniform interior ball condition implies the interior ball condition and that both 
the uniform interior and exterior ball conditions hold automatically for a C 2 -bounded domain. 

The Neumann boundary condition h is assumed to be continuous and uniformly bounded on dCl. 
Similarly, in the parabolic framework, the final-time data g is supposed to be continuous and uniformly 
bounded on Q. 



C 6 fe (fi) = \u e L°°(f7) : 30 D fi, domain, 3v € C^{0) s.t. u = v m 




It is equipped with the norm ||-|| C £(a) given by ||^)|| C fe(n) = X^- '^'^^)- 

i=0 



x = proj oa (x). 
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The real- valued function / in Ijl.ip is denned on K x !! x 1 x x S^. It is assumed throughout to 
be a continuous function of all its variables, and also that 

• / is monotone in T in the sense that 

(1.4) f{t,x,z,p,Ti+r 2 )<f(t,x,z,p,T 1 ) forr 2 >0. 

In the time-dependent setting (jl.ljl we permit / to grow linearly in \z\ (so solutions can grow exponentially, 
but cannot blow up). However we require uniform control in x (so solutions remain bounded as ||x|| — > oo 
with t fixed). In fact we assume that 

• / has at most linear growth in z near p = 0, T = 0, in the sense that for any K we have 

(1.5) \f(t,x,z,p,T)\<C K (l + \z\), 

for some constant Ck > 0, for all x £ and t, z £ M, when |j (j>, T)\\ < K. 

• / is locally Lipschitz in p and T in the sense that for any K we have 

(1.6) \f(t, x, z, P , r) - f(t, x, z, P ', r')| < c K (i + \z\)\\(p, r) - ( P ', r')|| , 

for some constant C K > 0, for all x £ O and t, z £ M, when || (p, T)\\ + \\ (p 1 , r')|| < K. 

• / has controlled growth with respect to p and T, in the sense that for some constants q, r > 1, 
C > 0, we have 

(1.7) |/(t,x,^,p,r)|<c(i + |z| + ||p||' + ||r|r r ), 

for all t, x, z,p and V . 



In the stationary setting fj 1 . 2|) our solutions will be uniformly bounded. To prove the existence of such 
solutions we need the discounting to be sufficiently large. We also need analogues of (|1.6[) - (|1.7[) but they 
can be local in z since z will ultimately be restricted to a compact set. In fact, we assume that 

• There exists r\ > such that for all K > 0, there exists C* K > satisfying 
(1.8) \f(x,z,p,T)\<(\-r))\z\+C* K , 

for all x £ fl, z € R, when < K; here A is the coefficient of u in the equation (jl.2|) . 

• / is locally Lipschitz in p and T in the sense that for any K and L we have 

(i-9) |/0r, z, P , r) - f(x, z, P ', r')| < c KfL \\(p, r) - ( P ', r')||, 

for some constant Ck,l > 0, for all x £ CI, when ||(p, r)|j + ||(p',r')|| < K and \z\ < L. 

• / has controlled growth with respect to p and T, in the sense that for some constants q, r > 1 
and for any L we have 

(i-io) |/( a; ,z,p,r)|<C7 L (i + l|p|| 9 + l|r|r), 

for some constant Cl > 0, for all x, p and T, and any |z| < L. 



2. The games 

This section present our games. We begin by dealing with the linear heat equation. Section [T"T] adresses 
the time-dependent problem depending non linearly on u; our main rigorous result for the time-dependent 
setting is stated here (Theorem 12.41) . Section 11.21 discusses the stationary setting and states our main 
rigorous result for that case (Theorem 12.71) . 
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2.1. The linear heat equation. This section offers a deterministic two-persons game approach to the 
linear heat equation in one space dimension. More precisely, let a < c and Q =]a,c[. We consider the 
linear heat equation on f2 with continuous final time data g and Neumann boundary condition h given 
by 

! u t + u xx — 0, for x E fl and t < T, 
^(x,t) = h(x), for xedn = {a,c} and t < T, 
on _ 
u(x, T) = g(x), for x £ f2 and t = T. 

Our goal is to capture, in the simplest possible setting, how a homogeneous Neumann condition can be 
retrieved through a repeated deterministic game. The game discussed here shares many features with 
the ones we will introduce in Sections 12. 2H2.3[ though it is not a special case. In particular, it allows to 
understand the way we need to modify the rules of the pioneering games proposed by Kohn and Serfaty 



in 2l| in order to model the Neumann boundary condition. 

There are two players, we call them Mark and Helen. A small parameter e > is fixed as are the 
final time T, "Helen's payoff" (a continuous function g: [a, c] — > R) and a "coupon profile" close to the 
boundary (a function h: {a, c} — >• R). The state of the game is described by its "spatial position" x G Q 
and "Helen's score" y E M. We suppose the game begins at time to- Since time steps are increments of 
e 2 , it is convenient to assume that T — to = Ke 2 , for some K. 

When the game begins, the position can have any value xo € f^; Helen's initial score is yo = 0. The 
rules are as follows: if, at time tj = to + je 2 , the position is Xj and Helen's score is yj, then 

• Helen chooses a real number pj. 

• After seeing Helen's choice, Mark chooses bj = ±1 which gives an intermediate position ij+i = 
xj + Aij where 

Axj = V2ebj € R. 

This position ij+i determines the next position Xj+i — Xj + Axj at time tj+% by the rule 

x j+ i = proj n (% + i) € tt, 

and Helen's score changes to 

(2.2) y j+1 = yj + Pj Axj ~ \\x j+1 ~ x j+1 \\h(xj + Axj ) . 

• The clock moves forward to tj+i = tj + e 2 and the process repeats, stopping when tx = T. 

• At the final time tx = T a bonus g{xn) is added to Helen's score, where xk is the final-time 
position. 

Remark 2.1. To give a sense to (|2.2p for all Axj, the function h, which is defined only on {a,c}, can be 
extended on ]a,c[ by any function Q — > R since \\xj+i — Xj+i\\ is different from zero if and only if ij+i ^ fi. 
Moreover, by comparing the two moves Axj and Axj, it is clear that \\xj+\ — = ||Ai£j — Aij||. 

Helen's goal is to maximize her final score, while Mark's goal is to obstruct her. We are interested in 
Helen's "value function" u E (xo,to), defined formally as her maximum worst-case final score starting from 
Xo at time to- It is determined by the dynamic programming principle 

(2.3) u £ (x, tj) — max min [u £ (x + Ax, tj+\) — pAx + \\Ax — Ax\\h(x + Ax)} , 

p£]R b—izl 

where Ax = \f2eb and Aa; = proj^(a; + Ax) — x, associated with the final-time condition 

u%x,T) =g{x). 

Evidently, if t = T - Ke 2 then 
(2.4) 

r 

u £ (xo,To) = max min • • • max min < g{xx) + —V2ebjPj + \\A£j — Axj\\h(xj + Axj) 
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where Axj — y/2ebj and Axj — pro^(xj + Axj ) — Xj . In calling this Helen's value function, we are using 
an established convention from the theory of discrete-time, two person games (see e.g. (l5|). 

By introducing the operator L e defined by 

(2.5) L e [x, 4>] = max mm [<j>(x + Ax) - pAx + \\Ax - Ax\\h(x + Ax)} , 

where Ax = \f2eb and Ax — proj^-(a; + Ax) — x, the dynamic programming principle (|2.3|) can be written 
in the form 

(2.6) u £ (x,t) = L E [x,u £ (-,t + e 2 )}. 

We now formally argue that u e should converge as e — > to the solution of the linear heat equation (|2.1[) . 
The procedure for formal passage from the dynamic programming principle to the associated PDE is 
familiar: we suppress the dependence of u e on e and we assume u is smooth enough to use the Taylor 
expansion. The first step leads to 

(2.7) u(x,t) w L e [x,u(-,t + e 2 )}. 

For the second step we need to compute L £ for a C 2 -function </>. By the Taylor expansion 

cj)(x + Ax) = 4>(x) + 4> x {x)Ax + -(j) xx (x)(Ax) 2 + 0(e 3 ) 

= 4>{x) + 4> x {x)Ax + \\Ax - Ax\\<p x (x)n(x) + ^<j) xx {x){Ax) 2 + 0(e 3 ), 

where x = projg^x), Ax — Ax = \\Ax — Ax\\n{x) with n defined on dQ by nix) = 1 if x = c and 

n(x) = — 1 if x = a. Substituting this expression in (|2.5[) . we deduce that for all C 2 -function 4>, 

(2.8) 

\o(s 2 ). 



L £ [x, <j>\ — 4>(x) + max min 



b=±i 



b x -p)Ax + ]^<f) xx (Ax) 2 + || Ax - Ax\\{h(x + Ax) - n{x)4> x ) 



It remains to compute the max min. If d[x) > s/2e, we always have Ax = Ax = V2eb, so that 
the boundary is never crossed and we retrieve the usual situation detailed in |2ll . Section 2.1]: Helen's 
optimal choice is p — (j) x and L E [x, <fi] — 4>{x) + £ 2 4> xx (x) + o(e 2 ). If d(x) < y2e, we still have Ax = y2be 
but there is a change: if the boundary is crossed, Ax — d(x) and || Ax — Ax\\ — \/2e — d(x). Suppose that 
Helen has chosen pel. Considering the min in (|2.8|) . Mark only has two possibilities b 6 {±1}- More 
precisely, suppose that x is close to c so that x — c and n(x) = 1; the case when x is close to a is strictly 
parallel. If Mark chooses 6=1, the associated value is 

V p ,+ = V2(<j) x - P )e + )^<t> xx d 2 {x) + (V2e - d(x))(h(c) - (f> x ), 
while if Mark chooses b = — 1, the associated value is 

V p - = -V2(<f> x ~ p)e + 4> xx e 2 . 
To determine his strategy, Mark compares V Pl - to V Pt +. He chooses b = — 1 if V Pl - < V Pt +, i.e. if 

V2(cj) x -p)e + ^(f> xx d 2 (x) + (V2e - d(x))(h(c) - <j) x ) > -V2(cp x - p)e + (j) xx e 2 , 
that we can rearrange into 

2V2(0 X - p)e > ^ xx (e 2 - - (V2e - d(x)) [h(c) - tf> g ] . 

This last inequality yields an explicit condition on the choice of p previously made by Helen 

1 I \ d(x)\ r , . . , . 1 / d 2 (x) 
(2-9) p < Popt := <t>* + r 1 - - V ^( c ) " W + Ts* 



2 V V2eJ L w J 2V2 V 2 £ 
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Meanwhile Mark chooses b = 1 if V Pl + < V Pl -, which leads to the reverse inequality p > p op t- The 
situation when V Pt + = V Pt - obviously corresponds to p — p op t- We deduce that 



L e [x,t 



max V p _ , V p . _ , max V p + 

P<Popt ' ' P>Popt 



Helen wants to optimize her choice of p. The functions V Pl + and V Pl - are both affine on 4> x — p. The first 
one is decreasing while the second is increasing with respect to p. As a result, we deduce that Helen's 
optimal choice is p — p opt as defined in (|2.9[) and L e [x,<fr] = Vp opt .+ = VJ, opt! _. We notice that Helen 
behaves optimally by becoming indifferent to Mark's choice; our games will not always conserve this 
feature, which was observed in [21|. Finally, for all C 2 -function </>, we have 

< 2 - 10 > 

[e 2 cb xx (x) +o(e 2 ), ifd(x)>V2e. 



o(e 2 ) and we formally derive the PDE by plugging (|2.f 01) in (|2.7p . This gives 

1 



Since u is supposed to be smooth, the Taylor expansion on t yields that u(-,t + e 2 ) = u(-, t) + ut(-,t)e 2 + 

in CHI). This 

(2.11) 0^ t+ < Vn^^)^"^ if^)<V2e, 
e 2 u xx +o(e 2 ), ifd(x)>A/2£. 

If x € fi, for e small enough, the second alternative in (|2.11[) is always valid so that we deduce from the 
£ 2 -order terms in (|2.1ip that u t + u xx = 0. If x is on the boundary dSl, then d(x) = 0, x = x and the first 
possibility in (|2.1ip is always satisfied. We observe that the £-order term is predominant since e e 2 . 
By dividing by £ and letting e — > 0, we obtain h(x) — u x {x) ■ n(x) = 0. 

Now we present a financial interpretation of this game. Helen plays the role of a hedger or an investor, 
while Mark represents the market. The position x is a stock price which evolves in as a function 
of time t, starting at xq at time to and the boundary dfl plays the role of barriers which additionally 
determine a coupon when the stock price crosses dtt. The small parameter e determines both the stock 
price increments Ax < V2~£ and the time step £ 2 . Helen's score keeps track of the profits and losses 
generated by her hedging activity. 

Helen's situation is as follows: she holds an option that will pay her g(x(T)) at time T (g could be 
negative). Her goal is to hedge this position by buying or selling the stock at each time increment. She 
can borrow and lend money without paying or collecting any interest, and can take any (long or short) 
stock position she desires. At each step, Helen chooses a real number pj (depending on Xj and tj), then 
adjusts her portfolio so it contains — pj units of stock (borrowing or lending to finance the transaction, so 
there is no change in her overall wealth). Mark sees Helen's choice. Taking it into account, he makes the 
stock go up or down (i.e. he chooses bj = ±1), trying to degrade her outcome. The stock price changes 
from Xj to xj + i = proj^Xj + Aij), and Helen's wealth changes by —\[2ebjPj + \\Axj — Axj\\h(xj + Axj) 
(she has a profit if it is positive, a loss if it is negative). The term ||AiEj — Axj\\h(xj + Axj) is a coupon 
that will be produced only if the special event Axj ^ fi happens. The hedger must take into account 
the possibility of this new event. The hedging parameter pj is modified close to the boundary but the 
hedger's value function is still independent from the variations of the market. At the final time Helen 
collects her option payoff g(xx)- If Helen and Mark both behave optimally at each stage, then we deduce 
by that 

K-l 

u £ (xo,to) + VZsbjPj — \\Axj — Axj\\h(xj + Axj) = g(xj<)- 
j=o 

Helen's decisions are in fact identical to those of an investor hedging an option with payoff g(x) and 
coupon h(x) if the underlying asset crosses the barrier dfl in a binomial-tree market with Ax = ^/2e at 
each timestep. 
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2.2. General parabolic equations. This section explains what to do when / depends on Du, D 2 u and 

also on u. We also permit dependence on x and t, so we are now discussing a fully-nonlinear (degenerate) 
parabolic equation of the form 

{d t u - f(t, x, u, Du, D 2 u) = 0, for x € and t < T, 
(Du(x, t), n(x)) = h(x), for x € dil and t < T, 

u(x,T) = g(x), for x 6 f2, 

where f2 is a C 2 -domain satisfying both the uniform interior and exterior ball conditions and the boundary 
condition h and the final-time data g are uniformly bounded, continuous, depending only on x. 

There are two players, Helen and Mark; a small parameter e is fixed. Since the PDE is to be solved 
in fi, Helen's final-time bonus g is now a function of x € f2 and Helen's coupon profile h is a function of 
.t € dfl. The state of the game is described by its spatial position x € f2 and Helen's debt z£l. Helen's 
goal is to minimize her final debt, while Mark's is to obstruct her. 

The rules of the game depend on three new parameters, a, /3, 7 > whose presence represents no loss 
of generality. Their role will be clear in a moment. The requirements 

(2.13) a<l/3, 
and 

(2.14) a + 0<l, 2a + 7 <2, max(/3 9 , fir) < 2, 

will be clear in the explanation of the game. However, the proof of convergence in Section [3] and consis- 
tency in Section U] needs more: there we will require 

(2.15) j<l-a, /3(q-l) <a + l, j(r-l)<2a, -yr < 1 + a. 

These conditions do not restrict the class of PDEs we consider, since for any q and r there exist a, f3 and 
7 with the desired properties. 

Using the language of our financial interpretation: 

a) First we consider U e (x, z,t), Helen's optimal wealth at time T, if initially at time t the stock 
price is x and her wealth is — z. 

b) Then we define u £ (x,t) or v £ (x,t) as, roughly speaking, the initial debt Helen should have at 
time t to break even at time T. 

The proper definition of U £ (x, z, t) involves a game similar to that of Section |2~T1 The rules are as follows: 
if at time tj = to + je 2 , the position is Xj and Helen's debt is Zj , then 

(1) Helen chooses a vector pj 6 R N and a matrix Tj € § N , restricted by 

(2.16) \\ Pj \\ <e-M|r,|| <e~\ 

(2) Taking Helen's choice into account, Mark chooses the stock price Xj+i so as to degrade Helen's 



outcome. Mark chooses an intermediate point ij+i = Xj + A Xj e R N such that 



(2.17) ||A%|| < e 1 -". 

This position %+i determines the new position Xj+i — Xj + Axj € f2 at time tj-\-\ by the rule 

(2.18) x j+1 = projjjOry+i). 

(3) Helen's debt changes to 

(2.19) Zj+i = Zj + pj ■ Axj + — (TjAxj, Axj) + e 2 f(tj,Xj , Zj,pj, Tj) — \\Axj — Axj\\h(xj + Axj). 

(4) The clock steps forward to tj+i = tj + e 2 and the process repeats, stopping when tx = T. At 
the final time Helen receives g(xn) from the option. 
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Figure 1. Rules of the game, near the boundary and inside the domain. 



This game is well-posed for all e > small enough. As mentioned in the introduction, the uniform 
exterior ball condition holds automatically for a C 2 -bounded domain. In this case, by compactness of 
dil, there exists e» > such that proj^ is well-defined for all x £ f2 such that d(x) < £*. It can be noticed 
that an unbounded C 2 -domain, even with bounded curvature, does not generally satisfy this condition. 
Since the domain Q satisfy the uniform exterior ball condition given by Definition 11.21 for a certain r, the 
projection is well-defined on the tubular neighborhood {x £ ~R N \Cl, d(x) < r/2} of the boundary. 

Remark 2.2. To give a sense to (|2.19p for all Axj, the function h which is defined only on the boundary 
can be extended on £1 by any function fl — > R since \\xj + i — is different from zero if and only if 

ij+i ^ f2. Moreover, by comparing Axj and Axj, one gets the relation 

Xj+\ = Xj+i + Axj — Aij. 

If £ Q, then Xj + \ = .ij+i and the rules of the usual game pi!/ are retrieved. Figure^ presents the 
two geometric situations for the choice for Mark: B(x,£ l ~ a ) C fl or not. 

Helen's goal is to maximize her worst-case score at time T, and Mark's is to work against her. Her 
value function is 

(2.20) U £ (xQ,z ,to) = max [gixx) - zk] ■ 

Helen's choices 

It is characterized by the dynamic programming principle 

(2.21) U e (x, z, tj) — maxmin U E (x + Ax, z + Az, ij+i) 

p,T Ax 

together with the final-time condition U e (x, z, T) = g(x) — z. Here Ax is Xj + i — Xj, Ax is determined by 

(2.22) Ax = Xj+i — Xj = projjy(:rj + Aij) — Xj, 

and Az = Zj+i — Zj is given by ()2.19j) . and the optimizations are constrained by (|2.16D and (|2.17|) . It is 
easy to see that the max/min is achieved and is a continuous function of x and z at each discrete time 
(the proof is by induction backward in time, like the argument sketched in [21]). 

When / depends on z, the function z h- > U e (x,z,t) can be nonmonotone, so we must distinguish 
between the minimal and maximal debt with which Helen breaks even at time T. Thus, following (ll| . 
we define 

(2.23) u £ (x , t ) = sup{z : U £ (x Q ,z a , t ) > 0} 
and 



(2.24) 



v £ (x , t ) = inf{z : U e (x ,z , t ) < 0}, 
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with the convention that the empty set has sup = — oo and inf = oo. Clearly v £ < u £ , and u £ (x,T) — 
v £ (x,T) = g{x). Since the definitions of u £ and v £ are implicit, these functions can not be characterized 
by a dynamic programming principle. However we still have two "dynamic programming inequalities". 

Proposition 2.3. Ifu £ (x,t) is finite then 

(2.25) u £ (x,t) < sup inf [u £ {x + Ax, t + e 2 ) 

- (^p ■ Ax + i (TAx, Ax) + e 2 f(t ) x,u s (x,t),p,T) - \\Ax - Ax\\h(x + Ax) 

Similarly, ijv £ {x,t) is finite then 

(2.26) v £ (x,t) > sup inf [v £ (x + Ax, t + e 2 ) 

p,T A * 

p ■ Ax + i (rAx, Ax) + e 2 f(t,x,v £ (x,t),p,T) - ||Ax - Ax\\h(x + Ax) 
The sup and inf are constrained by (|2.16p and (|2.17p and Ax is determined by (|2.22p . 

Proo f. The argument follows the same lines as the proof of the dynamic programming inequalities given 
in |2ll . Proposition 2.1]. R>r sake of completeness we give here the details. To prove (|2.25l) . consider 
z = u e {x,t). By the definition of u e (and remembering that U e is continuous) we have U e (x,z,t) = 0. 
Hence writing (|2.2ip . we have 

= maxminC/ 6 fx + Ax,z+p- Ax + - (TAx, Ax) + e 2 f{t,x, z,p,T) - ||Ax- Ax\\h(x + Ax),t + e' 
p,r Ax V 2 

We conclude that there exist p, T (constrained by (|2.16[0 such that for all Ax constrained by (I2.17[) . 
determining Aa; by (|2.22p . we have 

U £ (x + Ax,z+p - Ax + ^(rAx, Ax) + e 2 f(t,x, z,p,T) - ||Ax - Ax\\h(x + Ax),t + e 2 ^j > 0. 
By the definition of u e given by (|2.23p , this implies that 

z+p ■ Ax+ - (TAx, Ax) +e 2 f(t,x,z,p,T) - ||Ax- Ax\\h(x + Ax) <u £ (x + Ax,t + e 2 ). 
In other words, there exist p,T such that for every Ax, determining Ax by (|2.22p . 

z < u e {x + Ax,t + e 2 ) - (p ■ Ax + i(rAx, Ax) + e 2 f(t, x, z,p, T) - \\Ax - Ax\\h(x + Ax) 

Recalling that z — u e (x,t) and passing to the inf and sup, we get (|2.25l) . The proof of (I2.26P follows 
exactly the same lines. □ 

To define viscosity subsolutions and supersolutions, we shall follow the Barles and Berthame proce- 
dure [8], let us recall the upper and lower relaxed semi-limits defined for (t, x) G [0, T] x fi as 

(2.27) u(x,t) = limsup u e (y, tj) and v(x,t) = lira irtf_v £ (y,tj), 

tj-H f j — *■* 

where the discrete times are tj = T — je 2 . We shall show, under suitable hypotheses, that v and u are 
respectively viscosity super and subsolutions of p,12p . Before stating our rigorous result in Section [2. 2. 21 
the next section presents the heuristic derivation of the PDE (|2 . 1 2|) through the optimal strategies of 
Helen and Mark. 
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2.2.1. Heuristic derivation of the optimal player strategies. We now formally show that u e should converge 
as e — > to the solution of (|2.12| . Roughly speaking, the PDE (|2.12j) is the formal Hamilton Jacobi 
Bellman equation associated to the two-persons game presented at the beginning of the present section. 
The procedure for formal derivation from the dynamic programming principle to a corresponding PDE is 
classical: we assume u £ and v e coincide and are smooth to use Taylor expansion, suppress the dependence 
of if and v e on e and finally make e — > 0. That has already been done for x far from the boundary in [2]], 
Section 2.2] for / depending only on (Du, D 2 u). We now suppose that x is close enough of the boundary 
so that x can be nontrivial. By assuming u e = v e as announced and suppressing the dependence of u e 
on e, the two dynamic programming inequalities (|2.25p and (|2.26l) give the programming equality 

(2.28) u(x, t) k, sup inf [u(x + Ax, t + e 2 ) 

p,T A * 

- (p ■ Ax + i (TAx,Ax) + e 2 f(t,x,u(x,t),p,T) - \\Ax - Ax\\h(x + Ax) 

Remembering that Ax is small, if u is assumed to be smooth, we obtain 
u(x + Ax, t + e 2 ) + \\Ax - Ax\\h(x + Ax) 

« u(x,t) + e 2 u t + Du ■ Ax + ^ (D 2 uAx,Ax) + \\Ax - Ax\\h{x + Ax) 

« u(x, t) + e 2 u t + Du ■ Ax + || Ax - Ax\\ {h(x + Ax) - Du ■ n(x + Ax)} + ^ (D 2 uAx, Ax) , 

Ax - Ax 

since the outer normal can be expressed by nix + Ax) — : — - if x e O. Substituting this 

||Ax-A.t|| ^ 

computation in (|2.28p . and rearranging the terms, we get 

(2.29) w e 2 u t + maxmin [(Du - p) ■ Ax + \\Ax - Ax\\{h(x + Ax) - Du ■ n(x + Ax)} 

p,f Ax 

+ ^(D 2 uAx,Ax) - i (TAx,Ax)-e 2 f(t,x,u, P ,T) , 

where u, Du, D 2 u are evaluated at (x,t). We have ignored the upper bounds in (|2.16[) since they allow 
p, r to be arbitrarily large in the limit s — > (we shall of course be more careful in Section 0]). 

If the domain 51 does not satisfy the uniform interior ball condition, f2 can present an infinity number of 
"neck pitchings" of neck size arbitrarily small. To avoid this situation, the uniform interior ball condition 
is used to impose a strictly positive lower bound on these necks. If x is supposed to be extremely close 
to the C 2 -boundary and ||A£|| < e 1_Q , the boundary looks like a hyperplane orthogonal to the outer 
normal vector n(x), where x is the projection of x on the boundary dil (see Figure [2]). By Gram-Schmidt 
process, we can find some vectors e%, ■ ■ ■ ,e^ such that (ei = n(x), ei, ■ ■ ■ , e^r) form an orthonormal basis 
of WL N . In this basis, denote 





( 


• ' ' \^li)2<i<N ' 




(2.30) p = p!n(x)+p and V = {(Te,, ej)) 1<i]<N = 


{Fn)2<i<N 


f 






V : 







where pi e E, p G V 1 - = span(e 2 , • • • , ejv) and T = ((re», ej)) 2 <ij<jv e § 
Let us focus on the Neumann penalization term in (|2.29[) denoted by 



. . ., . . ,, . . . , , . , \ h(x + Ax) — Du(x) ■ nix + Ax), if x 4. f2, 

Fix) = Ai - Ax\\m(Ax) with mi Ax ) = I y ' v ; v h * _' 

V ' 11 11 V 1 V ' \m{Ax), if a GO, 

where m(Ax) is extended for x £ Vi by any function fh(Ax) (see Remark 12. 2|l . This contribution is 
favorable to Helen, P(x) > 0, if m(x) > 0, or to Mark, P(x) < 0, if m(x) < 0, and its size depends 
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|| Ax - Ax| 



Figure 2. Formal derivation for x near the boundary <90, notation: x — projg^x). 



on the magnitude of the vector Ax — Ax. Our formal derivation is local and essentially geometric, in 
the sense that our target is to determine the optimal choices for Helen by considering all the moves 
Ax that Mark can choose. By continuity of h and smoothness of u, the function m(Ax) is close to 
m = h(x) — Du{x) ■ n(x) if x ^ f2. We shall assume here that m(Ax), which serves to model the Neumann 
boundary condition, is locally constant on the boundary and equal to m. This hypothesis corresponds 
in the game to assume that in a small neighborhood, crossing the boundary is always favorable to one 
player. In order to focus only on the geometric aspects, this approach seems formally appropriate since 
it freezes the dependence of p(x) on m(x) by eliminating the difficulties linked to the local variations of 
m(x) like the change of sign. 

Hence, it is sufficient to examine 



(2.31) maxmin 

p,r Ax 



(Du - p) ■ Ax + || Ax - Ax||m + i (D 2 uAx, Ax) - ^ (TAx, Ax) - e 2 f (f, x, u,p, T) 



The formal proof will be performed in three steps. 

Step 1: To determine the optimal choice for Helen of p, we consider the e-order optimization problem 
M obtained from (j2.3ip by neglecting the second e-order terms 

(2.32) M = maxmin [(Du - p) ■ Ax + \\Ax - Ax\\m] . 

p Ax 

By writing Ax = (Ax)in(x) + Ax with Ax £ V 1 - and observing that || Ax — Ax|| depends only on (Ax)i, 
we decompose the max min (|2.32p into 



M = max min 

pi,p Ax 



(Du-p) ■ Ax + (Du! -pi)(Ax)i + || Ax - Ax|| 



max mm 

Pi KA^il^s 1 - 



(Du\ — pi)(Ax)i + || Ax — Ax||to + max min 

P HAxll^y^^ 



_(Du — p) ■ Ax 



-|(Ai)i| 2 



Noticing that the choices of p and p\ are independent from each other, we can successively solve the 
optimization problems. First of all, in order to choose p, let us determine 



M = max min 

p ||Ai||< v /e 2 - 2 <* 



_(Du — p) ■ Ax. 



-|(A£)i 



If Ax = ie 1 a n(x), Ax = and the min is always zero: Helen's choice is irrelevant. Otherwise, Helen 
should take p — projyi Du — Du, since otherwise Mark can make this max min strictly negative and 
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minimal by choosing Ax = —y/e 2 ~ 2a — |(Ax)i| 2 ^-p^%r with Ax ^ ±e 1 ~ a n(x). Thus Helen chooses 
p = Du, M = and M reduces to 

(2.33) M = maxmin [((Dm)i - p x )(Ax)i + ||Ax - Ax||m] . 

pi Ax 

To determine the remaining coordinate pi = p ■ n(x) of p, we now consider the optimization problem 
(12.33[) by restricting the possible choices made by Mark to the moves Ax which belong to the subspace 
V = Mn(x). Since || Ax|j < e x ~ a and Ax 6 V, we use the parametrization Ax = \e 1 ~ a n(x), A 6 [—1, 1]. 
If x G f2, the boundary is not crossed and ||Ax — Ax| = 0, while if x ^ f2 the boundary is crossed and 
|j Ax — Ax\\ = Ae 1- " — d(x). The intermediate point x — x £ dil separating the two regions corresponds to 
Ao = jr=L and || Ax — Ax\\ — 0. As a result, to compute the min in (|2.33l) . we shall distinguish these two 
regions by decomposing the global minimization problem into two minimization problems respectively on 
each region 

(2.34) M = max k(s p ) with n(s p ) = min(jV[i(s p ), M 2 (s p )), 
where s p = (Du — p) ■ n(x) and 

(2.35) Mi(s p ) = min Mi(Ai) with Mi(Ai) = {s p + m)£ 1_a Ai - d(x)m, 

Xq < Ai < 1 

(2.36) M 2 (s p ) = min M 2 (A 2 ) with M 2 (A 2 ) = s p e 1 - Q A 2 . 

— l<A2<Ao 

For fixed p, the functions defining M\ and M 2 are affine and can easily be minimized separately: 

• If s p + m > 0, JVti(sp) is attained for Ai = Ao and Mi(s p ) = d(x)s p . 

• If s p + m < 0, Mi(sp) is attained for Ai = 1 and Mi(s p ) = e 1 ~ a s p + (e 1_Q — d(x))m. 

• If Sp > 0, M 2 (s p ) is attained for A 2 = —1 and M 2 (sp) = — e 1 ~ a s p . 

• If s p < , M 2 (s p ) is attained for A 2 = Ao and M 2 (s p ) = d(x)s p . 

Geometrically, A e {— 1,1, Aq} corresponds to three particular moves: Ax = ±e 1 ~"n(x) and Ax = 
d(x)n(x). We are going to distinguish several cases to compute the max min according to the sign of s p 
and m. First of all, let us assume that m is positive. 

(CI) If Sp > then s p + m > and the optimal choices are (Ai, A 2 ) = (Ao, —1). It remains to minimize 
between (|2.35[) and (|2.36[) . Taking into account that d(x) < e 1 ~ a and s p > 0, we get by the 
definition of k(s p ) given by (I2.34p that n(s p ) — min{d(x)s p , — e 1 ~ a s p } = —e 1 ~ a s p . 

(C2) If — m < s p < then (Ai, A 2 ) = (Ao, Ao) and k(s p ) = Mi(s p ) = M 2 (s p ) = d(x)s p . 

(C3) If Sp < -m < then (Ai, A 2 ) = (1, A ) and Mi(s p ) = e 1 " Q s p + (e 1 "" - d{x))m and M 2 (s p ) = 
d(x)s p . By multiplying the inequality s p < —m < by (e 1_a — d(x)), we get 

n(sp) — min{e 1_Q s p + (e 1 "™ — c?(x))m, d(x)s p } = d(x)s p . 



By combining cases (CI) (C3) we conclude that if m > 0, 

!e 1_Q s p + (e 1_a — d(x))m, if s p < —m, 
d(x)sp, if — m < Sp < 0, 

-e l - a s pi if Sp > 0. 

The max of k is zero and reached at the unique value s p = Du ■ n(x) — p\ = 0. Since p — Du by the 
previous analysis, we conclude in (|2.30|) that if m > 0, Helen's optimal choice is p = Du. 

Let us now suppose that m is negative. 

(C4) If s p < then s p + m < and the optimal choices are (Ai, A 2 ) = (1, Ao). By the definition of 
k(s p ) given by (I2.34p . we obtain 

(2.37) k(s p ) = mmie 1 -" s p + (e 1 "" - d(x))m, d{x)s p } = e 1_0, 8 p + (e 1 "" - d(x))m. 
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(C5) If s p > — to > then (Ai, A 2 ) = (A , —1) and Mi(s p ) = d(x)s p and M 2 (s p ) = — e 1 a s p . By the 



definition of k(s p ) given by (|2.34p . we obtain k(s p ) — min {d(x)i 



•p > 



(C6) If < s p < -m, then (A 1; A 2 ) = (1, -1) and Mi(s p ) = £ 1 - a s p + (e 1 



-s 1 a s p . By the definition of k(s p ) given by (|2.34p , we obtain 



d(x))m and M2(s p ) 



k(s p ) 



{e^sp + ie 1 



d{x))m, 



_l-a , 



The target for Helen is to maximize this minimum with respect to s p . Both functions intervening 
in the minimum are afHne: the first one is afhne, strictly increasing and is equal to (e 1_Q — 
d{x))m < for s p — and to d(x)m > for s p = —to whereas the second function is linear and 



strictly decreasing and is equal to me 1 a < for 



As a result, there is a unique s* 



such that these two functions are equal and this value precisely realizes the max of n on [0, - 



d{x))m 



This 



Thus, the best that Helen can hope corresponds toe 1 a s* + (e 1 a 
gives 

s * = {Du-p).n{x) = -\(l-^± 
We immediately check that s* € [0, — ■rf] , which implies the condition s* + m < ^to < 0. Thus, 
max n(sp) = — (e 1_Q — d(x))m is greater than the minimum obtained in (|2.37|) . 



By combining cases (C4) (C6) we conclude that if to < 0, 



«(s p ) 



+ (e 1 



d{x))m, if s p < s*, 



if Sp > s* . 



The max of k is equal to k(s*) and reached for s p = Du ■ n(x) —pi= s* . 

Let us give an intermediate conclusion: if to > 0, Helen chooses p — Du whereas if to < 0, she chooses 



(2.38) 



p = Du 



m 
~2 



1 - 



d{x) 



n{x). 



Step 2: We are now going to take into account the second order terms in e in the optimization 
problem. If m > 0, once Helen has chosen p = Du, the optimization problem (|2.3ip reduces to computing 



(2.39) 



max mm 
r as 



Ax - Aac||m + - (D 2 uAx, Ax) - - (TAx, Ax) - e 2 f (t, x, u, Du, r) 



Mark is going to choose Ax ■ n(x) < 0, because otherwise the first e-order quantity \\Ax — Ax\\m will be 
favorable to Helen. Then considering Ax ■ n(x) < 0, we have Ai; = Ax and by symmetry of the quadratic 
form associated to (D 2 u — T), the optimization problem ()2.31[) reduces to 



(2.40) max min 

r Ai-n(i)<0 



((D 2 



T) Ai, Ai) - e 2 / (f, x, u, Du, V) 
"1 



= e max mm 
r Ai 



—e 



({D 2 u-T)Ax,Ax) - f(t, x, u, Du, r) 



Helen should choose Y < D 2 u, since otherwise Mark can drive e~ 2 ((D 2 u — Y)Ax, Ax) to — oo by a suitable 
choice of Ai. Thus, the min attainable by Mark is zero and is at least realized for the choice Ai = 0. 
Helen's maximization reduces to 

max [ut — f(t, x, u, Du, Y)]. 

T<D 2 u 

Since the PDE is parabolic, i.e. since / satisfies (|1.4[) . Helen's optimal choice is Y — D 2 u and (|2.29D 
reduces formally to u t — f(t, x, u, Du, D 2 u) = 0. 

If to < 0, Helen must now choose Y. In fact, we are going to see that the choice of p\ = p ■ n(x) 
obtained at (|2.38ll can be slightly improved by taking into account the additional terms containing D 2 u 
and r. Suppose Helen chooses p such that {p — Du)]yx = (notice that our first order computation 
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(|2.38p fulfills this condition) and Mark chooses a move Ax* realizing the minimum on Ax in (|2 . 31[) . We 
consider two cases depending on Ax*. 

Case a: if Ax* E V , we can restrain the minimization problem to the moves Ax which belong to 
V , Ax = Ax. Thus, the optimization problem (|2.31[) reduces to computing 

\e- 2 ({D 2 u - f) Ax, Ax) - / (t, x, u, p, T) , 



Mi. 



e max mm 

r Axev 1 - 



where F 



Helen should choose T < D 2 u, since otherwise Mark can drive e 2 ({D 2 u — T)Ax, Ax) 
to — oo by a suitable choice of Ax. By repeating the same argument of ellipticity of / already used for 
m > 0, Helen's optimal choice is T = D 2 u. 

Case b: if Ax* ^ V 1 -, there exists an unit vector v orthogonal to n(x) such that Ax* £ span(n(x), v). 
Thus, we restrain the minimization problem on Ax given by (|2.3ip to the moves Ax which belong to the 
disk D — span(n(x), v) PI B(s 1 ~ a ). This gives the optimization problem Md given by 
(2.41) 

M B = max min (Du - p)i(Ax)i + ||Ax- Ax||m+ - (D 2 uAx, Ax) (rAx, Ax) — e 2 f (t, x, u,p, T) 

px .1' AieD L 2 2 

by taking into account that p = Du and T = D 2 u. Neglecting — s 2 f(t, x, z,p, T), we want to compute 

max min NY s„, T, Ax) with 

s p ,r Ax 

(2.42) jV(s pj T, Ax) = s p (Ax)i + || Ax - Ax||m + i (D 2 uAx, Ax) - ^ (rAx, Ax) . 

Since || Ax|| < e 1_Q , we parametrize the disk D by Ax = \e 1 ~ a n(x) + [ie 1 ~ a v with A 2 + /j 2 < 1. Notice 
that the calculation of T\ v :— (Tn(x), v) for v orthogonal to n(x) implies the computation of (Tn(x))\ v ± . 
If Mark chooses Ax such that A > Aq for which the boundary is crossed, 



1 



(2.43) N(s p , T, Ax) = (s p + to)£ 1 " q A - d(x)m + -d 2 (x)(D 2 u) n - -A 2 e 



2„2-2a 



111 



whereas for Ax such that A < Ao for which the boundary is not crossed, 

1 



d(x) 



(D 2 u) lv - ATi„ e 



-2-2a 



(2.44) N(s p , T, Ax) = s p e^ a X + -A 2 ((C 2 U ) n - r u )e 



M(D 2 u) lv -Y lv )e 2 - 2a . 

For fixed A, Mark will always choose [i so that the last term is negative and maximal which leads to 

f sgn(|^ (D 2 u) lv - XT lv )VT= >?, if A < A < 1, 
fL \sgn(\((D 2 u) lv -T lv ))VT^X 2 ~, if-l<A<A . 

The min of N(s p , T, Ax) on /i depends only on A = Ax • n(x) and will be denoted below by 3\f(s p , T, A). 
By virtue of (|2.43[) and (|2.44[) it corresponds, for A > Ao, to 

N(s p ,T,A) = {s p +m)e 1 - a X~d{x)m+^d 2 {x)(D 2 u) 11 -^X 2 e 2 - 2a T 11 -Vl - A 2 ^^{D 2 u) lv -\T lv 
and, for A < Ao, to 



N( Sp ,r,A) 



a X + iA 2 ((^ 2 u ) u - r u )e 2 - 2 " - lAlVT - ~X?\(D 2 u) lv - T lv 



_2-2a 



The second order terms containing D 2 u and T being a little perturbation for e > small enough compared 
to (e 1_Q — d(x))m for d(x) <C £ 1_Q , it is sufficient to consider the case [(C6)1 which led to (|238)) and 
(Ai,A2) = (1,-1) corresponding to Ax = ±e 1 ~ a n(x). Therefore, we are going to compare the moves 
close to the optimal choices Ax = ±e 1 ~ a n(x) previously obtained by considering only the first terms 
in the Taylor expansion. More precisely, we may assume A ~ ±1, which leads to making the change of 
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variables Ai = 1 — pi, X 2 = — 1 + p 2 and take — > for i = 1,2. After some computations, we get a 
Taylor expansion in /?;, i = 1,2, in the form 

(2.45) N(s p , r, 1 - Pl ) = {s p + m)e 1 ~ a - d{x)m + ^d 2 (x)(D 2 u)n - ~£ 2 - 2c T n 



and 

(2.46) N(s p ,T,-l + p 2 )-- 



d ( x ) /n 2 \ r 



_2-2a 1 _l-a 

£ + pie 



(s p + m) + i 



1 11 



0(e 2 - 2 >? /2 ) 



(D 2 u) lv - T lv 



e 2 - 2a + p 2 e 1 ' a (s p - ((D 2 u) n - T 11 )e L - a ) + 0(e 



2-2a 3/2 
Pi 



First of all, we are now going to focus on the 0-order terms on the pi,p2 variables. Dropping the next 
terms corresponds to the two particular moves Ai = ±e 1_a n(a;) (p\ = p 2 = 0). For fixed T, since these 
terms containing T have the same contributions, we can omit the dependence of UST(- , T, 1) and N(-, T, —1) 
on r. Then, by repeating the same arguments already used, there exists a unique s 2 realizing the max of 
min(3\f(-, r, 1), ZNT(-, T, — 1)) for which both functions are equal. After some calculations, we find 



(2.47) 



S*o =-r I 1 



d(x) 



1 ( 1 a d 2 (x) 



(D 2 u) n - 



If 771 < 0, Helen will finally choose 



(2. 



Po P t{x) = Du + 



1 



d(x) 



1 f,- a d 2 {x) 



{D 2 u)i 



77.(3;). 



To complete the analysis on the 0-order terms on the p\ , p 2 variables, we are now going to see how Helen 
must choose T\\. By conserving only the 0-order terms, we obtain 



M D « -(e 1 a - d(x))m + max 

2 Tn 



\{e 2 - 2a + d 2 (x))(D 2 u) n - ^ 2_2< T U - s 2 f (t, x, u,p opt ,T) 



Since Tn cannot counterbalance the first order term, the Tn-term and the second order terms are 
gathered. Helen wants to make the best choice, so she is going to choose Tn such that 



1 



2-2q 



1 



■d 2 (x))(D"u) u - ^-^Tn > 0. 



By ellipticity of /, Helen will choose fn such that this upper bound on Tn is attained. She takes 



(2.49) 



1 



d 2 (x) 



(D 2 u) 1 



,2-2q 



{s* 2 + m) Pl e i - a + 0(e 2 - 2a Pl ), 



It remains to determine Ti v . By plugging the optimal choices s 2 , corresponding to p pt, and Tn, respec- 
tively given by (|2T47j) and (|2T4^)l in (j2T4"5]) - (j2T4T)]) , we have 

N(^,r, 1 - p x ) = ^(e 1 -* - d{x))m - ^\^L{D 2 u)i v - T lv 
1 



N(«S,r, -1 + pi) = -{e l - a - d{x))m - ^2p 2 \{D 2 u) lv - T lv \ e 2 ~ 2a + s* 2 p 2 e L - a + 0(e 2 ~ 2a p 2 ). 

Dropping the 0(e 2 ~ 2a pi) terms and noticing that > and — (sj + m) > for e small enough, the two 

minimization problems minN^, T, 1 — pi), i G {1, 2} for Mark reduce to find 

Pi 

min f(p), where f(p) = a^fp + bp, 

0<p<l 



with a < < b. Differentiating /, the minimum of / is attained at \fff 



2b 



We can notice that this 



computation is equivalent to formally differentiating the Taylor expansion of (|2.43|) — (|2.44p . Conserving 
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the predominant terms and dropping the next terms, the minimum of !N(s 21 T, 1— pi) and INffsJ , T, 
are respectively attained at 



-P2) 



1 \^(D^ u ) lv -T lv \ 
V2 \s* 2 + m\ 



and 



1 \{D 2 u) lv -T lv \ 
V2 141 



Assuming formally that these approximations are in fact equalities, we obtain 



(2.50) 



(2.51) 



Pl) = \{e 1 - a 



d{x))m- i 



d{x) 2 



-3 — 3a 



+ 0(e 4 - 4a ), 



d(x))m - - \ (D 2 u)i v 



2 £ 



3-3a 



+ 0(e 4 - 4 "). 



Helen now has to choose r lt) such that min{jV(s2, I\ 1 — P*)> N(s2> I\ — 1 + p^)} is maximal. We could 
compute the optimal value of Ti v on the e 3_3a -terms. However, it is not very useful. Since to is a constant 
and e 3 ~ 3a <C e 2 by (|2.13p . the £ 3 ~ 3a -terms are negligible compared to —e 2 f(t, x, u,p opt , T) that we have 
omitted until now. For instance Helen can fix Ti^ such that one of the two terms depending on in 
(|2.50| and (|2.51| is equal to zero: = {D 2 u)\ v or T\ v — ^l-{D 2 u)i v . The two choices are equivalent 
because Mark can reverse his move Ai. For sake of simplicity, we assume Helen chooses T\ v = (D 2 u)i v . 
It is worth noticing that this expansion holds if m is far from zero and we shall modify our arguments 
very carefully in Section U when to is negative but small with respect to a certain power of e. 

Thus, if to < 0, Helen will choose 



(2.52) 



r opt (x) = d 2 u + 



-1 



d 2 (x) 

c 2-2a 



(D 2 u), 



Unlike the usual game [211 ], when Helen chooses p and T optimally, she does not become indifferent to 
Mark's choice of Ax. More precisely, it depends on the projection of Ax with respect to n(x). Our games 
always have this feature. 

Step 3: Now let us go back to the original optimization problem (|2.29p . If to = 0, by letting e 
we get h(x) — Du(x) ■ n(x) = 0. Otherwise, (|2.29t formally reduces to 
(2.53) 

-(e 1_Q - d(x))m - e 2 f(t,x,u,p op t{x),T opt (x)) 



0, 







e 2 u t 



—e 2 f(t, x, u, Du, D 2 u), 



-o(e 2 ), if d(x) < e 1 ' and to < 0, 
if d{x) > e 1 "" or m > 0, 

with p op t and r opt respectively defined by ()2.48|) and ()2.52p . If X G ft, for £ small enough, the 
second relation in (|2.53|) is always valid so that we deduce from the e 2 -order terms in ()2.53|) that 
u t — f(t,x,u,Du,D 2 u) = 0. If x £ <9f2, d{x) — and we distinguish the cases m > and to < 0. 
If to > 0, one more time the second relation in (|2.53[) is always valid so that u t — f(t, x, u, Du, D 2 u) = 0. 
Otherwise, if to < 0, the first relation in (|2.53p is always satisfied. We observe that the e-order term 



is predominant since e a 3> e . By dividing by e a and letting 



0. we obtain to = that leads 



to a contradiction since we assumed to < 0. Therefore, we have formally shown that on the boundary 
h(x) — Du(x) ■ n{x) — or u t — f(t, x, u, Du, D 2 u) = 0. 



2.2.2. Main parabolic result. We shall show, under suitable hypotheses, that u and v are respectively 
viscosity sub and supersolutions. A natural question is to compare u and v. This is a global question, 
which we can answer only if the PDE has a comparison principle. Such a principle asserts that if u is 
a subsolution and v is a supersolution then u < v. If the PDE has such a principle then it follows that 
u < v. The opposite inequality is immediate from the definitions, so it follows that u = v, and we get a 
viscosity solution of the PDE. It is in fact the unique viscosity solution, since the comparison principle 
implies uniqueness. 



Theorem 2.4. Consider the final-value problem (|2.12p where f satisfies (|1.4p - (|1.7p . g and h are con- 
tinuous, uniformly bounded, and ft is a C 2 -domain satisfying both the uniform interior and exterior ball 
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conditions. Assume the parameters a, j3, 7 satisfy (|2.13[) - (|2.15[) . Thenu and v are uniformly bounded on 
Q x [i, , T] for any < T, and they are respectively a viscosity subsolution and a viscosity supersolution 
of (|2.12[) . If the PDE has a comparison principle (for uniformly bounded solutions), then it follows that 
u £ and v £ converge locally uniformly to the unique viscosity solution of ([2.1211 . 

This theorem is an immediate consequence of Propositions 13.41 and 15.11 

In this theorem, we require the domain VL to be C 2 . This assumption is crucial for the proof of 



Proposition 13.41 case (iii) corresponding to the convergence at the final-time in the viscosity sense (see 
Remark l3.5p . It can also be noticed that it is this part of Proposition l3. 41 which allows to use a comparison 
principle for the parabolic PDE. On the other hand, since the game already requires the uniform interior 
and exterior ball conditions, the domain fl is in fact at least C ' . It remains an open question to overcome 
the analysis in this case. 

As mentioned in [2l| . some sufficient conditions for the PDE to have a comparison result can be found 
in Section 4.3 of (llj . In our framework, we can emphasise on the comparison principle obtained by Sato 



271 . Theorem 2.1] for a fully nonlinear parabolic equation with a homogeneous condition. The reader is 
also referred to the introduction for other references about comparison and existence results. Note that 
most comparison results require f(t,x,z,p,T) to be nondecreasing in z. 

We close this section with the observation that if U e (x, z,t) is a strictly decreasing function of z then 
v £ (x,t) — u e (x,t). A sufficient condition for this to hold is that / be nondecreasing in z: 

Lemma 2.5. Suppose f is non- decreasing in z in the sense that 

f(t, x, Zi,p, r) > f(t, x, zq,p, r) whenever Z\ > zq. 

Then U £ satisfies 

U £ (x, z\, tj) < U £ (x, zq, tj) — (zi — zq) whenever Z\ > zq, 
at each discrete time tj~T — je 2 . In particular, U £ is strictly decreasing in z and v £ = u £ . 

Proof. The whole space case is provided in |2ll . Lemma 2.4]. For our game, it suffices to add — ||Ai — 
Ax\\h(x + Ax) in the expressions of Azo and Azi defined in the proof of |2ll . Lemma 2.4]. The rest of 
the proof remains unchanged. □ 

2.3. Nonlinear elliptic equations. This section explains how our game can be used to solve stationary 
problems with Neumann boundary conditions. The framework is similar to the parabolic case, but one 
new issue arises: we must introduce discounting as in [2l|, to be sure Helen's value function is finite. 
Therefore we focus on 



(2.54) 



f(x,u,Du,D 2 u) + Xu = 0, in fl, 
(Du, n) — h, on dfl, 



where Q is a domain with C 2 -boundary and satisfies both the uniform interior and exterior ball condition 
presented in the introduction. The constant A (which plays the role of an interest rate) must be positive, 
and large enough so that (|1.8|) holds. Notice that if / is independent of z then any A will do. 

We now present the game. The main difference with Section |2~21 is the presence of discounting. The 
boundary condition h is assumed to be a bounded continuous function on dfl. Besides the parameters 
a, /3, 7 introduced previously, in the stationary case we need two new parameters, m and M, and a 



C 2 (i7)-function ip such that 



(2.55) 7^ = IWI«> + 1 ondn. 

on 

It suffices to construct ipx such that it is C 2 (fl) and satisfies = 1 on the boundary. Then we can 
define tp by ip = (||/i|joo + 1)^1- The existence and construction of such a function ipi for a C 2 -domain Q 
satisfying the uniform interior ball condition is discussed at the end of this section. 
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From m and tp we construct a function \ denned by 

(2.56) x{x) = m+ ||VHz,~(n) +ip(x). 

Both m and M are positive constants, which also yield that \ 1S positive. M serves to cap the score, 
and the function \ determines what happens when the cap is reached. We shall in due course choose m 
such that m + 2||V>||z°° = M — 1 and require that M is sufficiently large. Like the choices of a, (3, 7, the 
parameters M, m and the function ifj are used to define the game but they do not influence the resulting 
PDE. As in Section |2~21 we proceed in two steps: 

• First we introduce U E (x, z), the optimal worst-case present value of Helen's wealth if the initial 
stock price is x and her initial wealth is — z. 

• Then we define u E (x) and v E (x) as the maximal and minimal initial debt Helen should have at 
time t to break even upon exit. 

The definition of U E (x, z) for x € Q involves a game similar to that of the last section: 

(1) Initially, at time to = 0, the stock price is xq = x and Helen's debt is zq = z. 

(2) Suppose, at time tj = je 2 , the stock price is Xj and Helen's debt is Zj with \zj\ < M. Then Helen 
chooses a vector pj £ M. N and a matrix Tj G § N , restricted in magnitude by (|2 . 16|) . Knowing 
these choices, Mark determines the next stock price Xj+i = Xj + Ax so as to degrade Helen's 
outcome. The increment Ax allows to model the reflection exactly as in the previous subsections. 
Mark chooses an intermediate point ij+i = Xj + Aobj € R N such that 

||A4j|| < e 1 - 01 . 

This position Xj+i determines the new position + Axj at time tj+\ by 

Xj+i = proj TT (% +1 ). 
Helen experiences a loss at time tj of 

(2.57) Sj = pj ■ Aij + — (TjAxj, A£j) + e 2 f(xj, Zj,pj, Tj) — ||A^- — Axj\\h(xj + Axj). 
As a consequence, her time tj+x = tj + e 2 debt becomes 

Zj+i = e Xe \zj+5 3 ), 

where the factor e takes into account her interest payments. 

(3) If Zj+i > M, then the game terminates, and Helen pays a "termination-by-large-debt penalty" 
worth e Ae (x( x j)~ $j) a t time tj+\. Similarly, if Zj+i < —M, the the game terminates, and Helen 
receives a "termination-by-large-wealth bonus" worth e Ae2 (x( x j) + Sj) at time tj+\. If the game 
stops this way we call tj+i the "ending index" tx- 

(4) If the game has not terminated then Helen and Mark repeat this procedure at time tj +1 = tj +e 2 . 
If the game never stops the "ending index" tx is +00. 

Helen's goal is a bit different from before, due to the presence of discounting: she seeks to maximize the 
minimum present value of her future income, using the discount factor of e~° for income received at 
time tj. If the game ends by capping at time tx with Zk > M, then the present value of her income is 

U s (x ,z ) = -zo -S - e-^&i e-^-^dK-i - e^ K -^ (xfcjc-i) - 5 K -i) 

= e -{K-l)Xs> { _ ZK _ i _ x{xK _ i)) _ 

Similarly, if the game ends by capping at time tx with zx < —M, then the present value of her income 
is 

U%x , z ) = -zo -5 - e^^Sx-i + e^-^ixixx-i) + S K -i) 

= e- (K - 1 ^ 2 (-z K - 1 +x(xx-i)). 
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If the game never ends (since Z j and X (xj) are uniformly bounded), we can take K = oo in the preceding 
formula to see that the present value of her income is 0. 

To get a dynamic programming characterization of U £ , we observe that if |zo| < M then 

(e- Xs2 U s (x 1 ,z 1 ), i£\ Zl \<M, 
U £ (x , zq) = supinf < -z - xi^o), if zi > M, 
{-zo + x{®o), ifzi<-M. 

Since the game is stationary (nothing distinguishes time 0), the associated dynamic programming principle 
is that for \z\ < M, 



e 



U £ {x,z) 



U £ (x',z'), if \z'\<M, 

(2.58) U £ (x,z) = supinf { - z - x (x), if z' > M, 

n p Ax I 

1 X {x), if z 1 < -M, 

where x' — pro]^(x + Ax) and z' = e Xe (z + 5), with 5 defined as in (|2.57p . Here p, T and Aa; are 
constrained as usual by (|2.16[) - (|2.17j) . and we write sup / inf rather than max / min since it is no longer 
clear that the optima are achieved (since the right-hand side is now a discontinuous function of p, T and 
Ax). The preceding discussion defines U e only for \z\ < M; it is natural to extend the definition to all z 
by 

-z — x( x )i f° r z ^ M, 

-Z + Xi x ); f° r z < —M, 

which corresponds to play being "capped" immediately. Notice that when extended this way, U £ is strictly 
negative for z > M and strictly positive for z < —M. 

The definitions of u s and v £ are slightly different from those in Section 12.21 

(2.59) u £ (x ) = sup{2 : U £ (x ,z ) > 0}, 

(2.60) v £ (x ) = inf{z : U £ (x Q , z ) < 0}. 

The change from Section |2~21 is that the inequalities in (|2.23p - (|2.24p are strict. 

Proposition 2.6. Let mi, M be two constants such that < mi < M. Then whenever x € f2 and 
—mi < u £ (x) < M we have 

(2.61) u £ (x) < supinf \e- X£ \ £ {x + Ax) 

p ■ Ax + i {TAx, Ax) + e 2 f(x,u £ (x),p,T) - \\Ax - Ax\\h(x + Ax) 

for e small enough (depending on m\ and the parameters of the game but not on x). Similarly, if x G £1 
and —M < v £ (x) < mi then for e small enough 

(2.62) v £ (x) > supinf \e~ X£2 v £ (x + Ax) 

P ,r A * - 

p ■ Ax + i (TAx, Ax) + e 2 f(x,v £ (x),p,T) - \\Ax - Ax\\h(x + Ax] 

As usual, the sup and inf are constrained by (|2.16p and (|2 . 1 7|) and Ax is determined by (|2.22[) . 

Proof. We shall focus on (|2.61[) ; the proof for (|2.62l) follows exactly the same lines. Since —mi < u £ (x) < 
M, there is a sequence z k f u £ (x) such that U £ (x,z k ) > 0. Since u £ (x) is bounded away from —M, 
we may suppose that z k also remains bounded away from —M. Dropping the index k for simplicity of 
notation, consider any such z = z k . The fact that U £ (x, z) > tells us that the right-hand side of the 
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dynamic programming principle (|2.58l) is positive. So there exist p, T constrained by (|2.16l) such that for 
any Ax satisfying (|2.17|) . 

(e- Xe2 U s (x',z'), if \z'\<M, 

(2.63) < sup inf \ -z - x(x), if z' > M, 

v r Ax I 

{-z + x(x), itz'<-M, 

where x' = proj^(x + Ax) and z' = e Xe (z + 8). Capping above, the alternative z' > M , cannot happen, 
since otherwise we compute 

-z - x(x) = -e~ Xe ' V - 8 - x(x) < -Me~ Xe * - 8 - m < -8 - m < 0, 

for e small enough because 8 is bounded by a positive power of s. This sign is a contradiction to our 
assumption (|2.63|) . If e is sufficiently small, capping below (the alternative z' < —M) cannot occur either, 
because z is bounded away from — M and 8 is bounded by a positive power of e. Therefore only the first 
case can take place 

< U e {x + Ax,e Xe \z + 8)), 
whence by the definition of u e given by (|2.59p , we deduce that 

u e {x + Ax) > e Xe \z + 8). 
Thus, we have shown the existence of p, T such that for every Ax, 

(2.64) z < e~ A£ V(x + Ax) - ( p ■ Ax + -{TAx , Ax) + e 2 f(x 7 z 7 p, T) -\\Ax- Ax\\h(x + Ax) 



Recalling that z = z k t u e (x), we pass to the limit on both sides of (|2.64p . with p, T held fixed, to see 
that 

u e {x) < e~ Xs2 u £ (x + Ax) - ^p- Ax+ i(rAx, Ax) + e 2 /(x, u 6 {x),p, T) - \\Ax - Ax\\h(x + Ax) 

Since this is true for some p, T and for every Ax, we have established (|2.6ip . □ 

The PDE (|2.54|) is the formal Hamilton- Jacobi-Bellman equation associated with the dynamic pro- 
gramming inequalities (|2.61j) - (|2.62j) . by the usual Taylor expansion, if one accepts —M < v e « u £ < M. 
Rather than giving that heuristic argument which is quite similar to the one proposed in the parabolic 
setting, we now state our main result in the stationary setting, which follows from the results in Sections 
H]and[S] It concerns the upper and lower relaxed semi-limits, defined for any x G O, by 

(2.65) u{x) — limsupii e (y) and v (x) = lim inf v e (y) , 

y — yx V~* x 

with the convention that y approaches x from f2 (since u e and v 6 are defined on fl). 

Theorem 2.7. Consider the stationary boundary value problem (|2.54p where f satisfies (|1.4p and (|1.8p - 

(jl.lOp . g and h are continuous, uniformly bounded, and Q, is a C 2 -domain satisfying both the uniform 
interior and exterior ball conditions. Assume the parameters of the game a, ft, 7 fulfill (|2.13p - (l2.15p . ip € 
C 2 (il) satisfies (|2. 551) . x € C 2 (£!) is defined by (|2.56p . M is sufficiently large andm = M —l — 2\\ip\\ Loo( j^y 
Then u s and v E are well-defined when e is sufficiently small, and they satisfy \u £ \ < x an d \ v& \ 2^ X- Their 
relaxed semi-limits u and v are respectively a viscosity subsolution and a viscosity supersolution of (|2.54p . 
// in addition we have v < u and the PDE has a comparison principle, then it follows that u e and v e 
converge locally uniformly in f2 to the unique viscosity solution of (|2.54p . 

This is an immediate consequence of Propositions 13.61 and 15.51 A sufficient condition for v < u is that 
/ is nondecreasing in z. As mentioned in [21] , sufficient conditions for the PDE to have a comparison 
principle can be found for example in Section 5 of 12 1, and (for more results) in 
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Let us now go back to the existence and the construction of ip\ g C- 2 (il) such that 



fi = lon an, 

~ • - on 

that we will need at various points of the paper. If fl is of class C 2 and satisfies the uniform interior ball 
condition of Definition 11.21 for a certain r, d is C 2 on f2(3r/4) and an explicit suitable function is 



(2.66) 



ipi (x) 




d(x) 



if d(x) < r/2, 
if d(x) > r/2. 



It is clear that supp ipi C Q(r/2), ipi(ty c [0, 1] and is C 2 on f2(r/2). Then, for all x such that 
g?(x) = |, Dipt and D 2 V>i are continuous at x. Thus -0i is C 2 on f2. It is easy to check that the two first 
derivatives of ipi are also bounded and that = 1 on the boundary. Hence, the function ipi defined by 
(|2.66p has all the desired properties. 

Remark 2.8. If ft is a domain with C 2,a -boundary where a > 0, the Schauder theory fldi . Theorem 6.31] 
ensures the solution ip of the elliptic problem 

in f2, 
on dfl, 




is C 2,a (fl). In addition, the estimate || V'llc'^ffT) 

< C n (l 

depending only on the domain. 



\h\\L°°) holds for a certain constant Cf, 



3. Convergence 



This section presents our main convergence results .Jinking the limiting behavior of v e and u e as e — » 
to the PDE. The argument uses the framework of [10j and is basically a special case of the theorem 
proved there. 

Convergence is a local issue: in the time-dependent setting, Proposition 13.41 shows that in any region 
where the lower and upper semi-relaxed limits v and u are finite they are in fact viscosity super and 
subsolutions respectively. The analogous statement for the stationary case is more subtle. In fact, we 
will need global hypotheses on / at Section HT21 to ensure that u E and v £ are well-defined and satisfy the 
dynamic programming inequalities (|2.6ip - (|2.62p . Thus, we cannot discuss about v or u without global 
assumptions on /. 

3.1. Viscosity solutions with Neumann condition. Our PDEs can be degenerate parabolic, degen- 
erate elliptic, or even first order equations. Therefore, we cannot expect a classical solution, and we 
cannot always impose boundary data in the classical sense. The theory of viscosity solutions provides 
the proper framework for handling these issues. We review the basic definitions in our setting for the 
reader's convenience. We refer to [J], [l2| and [l6| for further details about the general theory. Consider 
first the final- value problem (|2.12l) in f2, 

—u t + f(t, x, u, Du, D 2 u) =0, for x e Vt and t <T, 
< (Du(x, t), n(x)) = h(x), for x e dfl and t < T, 

u{x,T) — g[x), for x G ft. 

where f(t,x,z,p,T) is continuous in all its variables and satisfies the monotonicity condition (|1.4I) in its 
last variable. We must be careful to impose the boundary condition in the viscosity sense. 

Definition 3.1. A real-valued lower- semicontinuous function u(x,t) defined for x € £1 and t„ < t < T is 
a viscosity supersolution of the final-value problem (|2.12| if 

(PI) for any (xo,to) with xq € O and t* < to < T and any smooth (f>(x,t) such that u — (f> has a local 
minimum at (xo,to), we have 

d t (f>(x ,to) - f(t , x , u(x , t ), D(f>(x , t ), D 2 4>(xo, to)) < 0, 
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(P2) for any (xo,to) with xq £ dfl and t* < to < T and any smooth <fi(x,t) such that u — <p has a local 
minimum on f2 at (xo,to), we have 

max{-(dt<f>(x ,to)-f(t ,xo, u(x , to), D<f>(x , t ), D 2 4>(x , to))), (D<f>(x , t ), n(x )) - h(x )} > 0, 
(PS) u > g at the final time t = T. 

Similarly, a real-valued upper- semicontinuous function u(x, t) defined for x £ f2 and f* < t < T is a 
viscosity subsolution of the final-value problem (|2.12p if 

(PI) for any (xo,to) with xq £ f2 and t* < to < T and any smooth <j)(x,t) such that u — <fi has a local 
maximum at (xo,to), we have 

d t 4>(x ,t ) - /(to, xq, u(x , t ), D(f>(x , to), D 2 (f>(x , t )) > 0, 

(P2) for any (xo,to) with xq £ dQ and t* < to < T and any smooth <f>(x,t) such that u — <fi has a local 
maximum on SI at (xo,to), we have 

min{-(<9 t </>(xo, t )-f(to,x , u(x , to), Dcj)(xo, to), D 2 cj)(xo, to))), (D(/>(xo, to), n(x )) - h(x )} < 0, 
(PS) u < g at the final time t = T . 

A viscosity solution of (|2.12|l is a continuous function u that is both a viscosity subsolution and a viscosity 
supersolution of ()2.12j) . 

In the stationary problem (|2.54[) . the definitions are similar to the time-dependent setting. 

Definition 3.2. A real-valued lower- semicontinuous function u(x) defined on fl is a viscosity supersolu- 
tion of the stationary problem (|2.54[) if 

(El) for any io 6 ^ and any smooth <f>(x) such that u — <fr has a local minimum at xq, we have 

f(x ,u(x ), D(f>(x ), D 2 (j)(x )) + Xu(x a ) > 0, 
(E2) for any xo € dQ and any smooth (j)(x) such that u— <f> has a local minimum on O at Xo, we have 
max{/(a;o, u(xo), D<f>(xo), D 2 4>(x )) + Xu(x ), (D<j>(xo),n(xo)) - h(x )} > 0. 

Similarly, a real-valued upper-semicontinuous function u(x) defined on CI is a viscosity subsolution of the 
stationary problem (|2.54| if 

(El) for any xq € and any smooth <p(x) such that u — <f> has a local maximum at xq, we have 

f(x ,u(x ), D(j)(xo), D 2 (j)(x )) + Xu(x a ) < 0, 

(E2) for any xq £ dfl and any smooth (f>(x) such that u — 4> has a local maximum on f2 at Xq, we have 

min{/(x , u(xq), D<f>(xo), D 2 (/)(xo)) + Xu(x Q ), (D<f>(xo), n(x )) - h(x )} < 0. 

A viscosity solution of (|2.54p is a continuous function u that is both a viscosity subsolution and a viscosity 
supersolution of (|2.54p . 

In stating these definitions, we have assumed that the final-time data g and the boundary Neumann 
condition h are continuous. In Definition 13. 11 the pointwise inequality in part |(P3"J| can be replaced by an 
apparently weaker condition analogous to part |(P2)| The equivalence of such a definition with the one 
stated above is standard, the argument uses barriers of the form <f>(x, t) = \\x — xq\\ 2 /S+ (T — t)/fi + Kd(x) 
with S and fx sufficiently small, and is contained in our proof of Proposition l3.4((rii) We shall be focusing 
on the lower and upper semi-relaxed limits of v e and u e , defined by (I2.27[) in the time-dependent setting 
and (|2.65[) in the stationary case. 

We now provide a key definition to deal with the Neumann boundary condition within viscosity so- 
lutions framework which will be essential all along the article. We introduce some applications which 
give bounds on the Neumann penalization term for a smooth function and x close to the boundary. This 
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approach is well-suited to the viscosity solutions framework. More precisely, we define the applications 
m e and M e , for all x G 0(e 1_Q ) and (ft G C 1 ^), by 

(3.1) m%[(j>]:= inf {h(x + Ax) - D(ft(x) ■ nix + Ax)} , 

x+Ax^Q, 
Ax 

(3.2) M*[(j)}:= sup {h(x + Ax) - D(ft[x) ■ n{x + Ax)} , 

x+A££Cl 
Ax 

where Ax is constrained by (|2.17p and determines Ax by (|2.22[) . Notice that the functions m e [0] and 
Af £ [0] are bounded by + ll-O^llL^fn)- Since h is supposed to be continuous, the following property 

clearly holds. 

Lemma 3.3. Let x G 90 and (ft G C 1 (0). Suppose there exists a sequence (sk,%k)kEN i n ^+ x 
convergent to (0,x) such that for all k large enough, Xk G f2(e^,~ Q ). Then 

lim m Xk \<p] = lim M Xk \<p] = h(x) - D(j>(x) ■ nix). 

fc->+oc fc fc->+oo * 

Similarly, let (ft G C 1 (0 x [0,T]). Suppose there exists a sequence (ek,Xk,tk)keti ^ n x O x [0,T] 
convergent to (Q,x,t) such that for all k large enough, Xk G r2(e^, _Q ). TTien 

lim m Xk M-> **)] = lim M X° [^0. *fc)l = ~ *) ' »(*)• 

3.2. The parabolic case. We are ready to state our main convergence result in the time-dependent 
setting. At first sight, the proof seems to use the monotonicity condition (|1.4p . The proof relies on 
the consistency of the numerical scheme, Propositions I4.5| 14.101 and I4.16| which are proved in Section 
|U Proposition 14.161 is necessary to deal with the degeneration of the consistency estimates due to the 
Neumann boundary condition. So we also require that f(t,x,z,p,T) satisfy (jl.6)l - (|1.7)l . and that the 
parameters a, (3, 7 satisfy (|2.13p - (|2.15p . 

Proposition 3.4. Suppose f and a, j3 } 7 satisfy the hypotheses just listed. Assume furthermore that u 
and v are finite for all x near Xq and all t < T near to . Then 

i. if to < T and xq G fl, then u is a viscosity subsolution at Xq and v is a supersolution at xq (i.e. 

each one satisfies part \(P1) \ of the relevant half of Definition \3.1\ at xo )■ 
ii. if to <T and xo G dfl, then u is a viscosity subsolution at xq and v is a supersolution at xo (i.e. 

each one satisfies part \(P2) \ of the relevant half of Definition \3.1\ at xq ). 
Hi. if to — T, then u(xq,T) = g(xo) and v(xo,T) = g(xo) (in particular, each one satisfies part \(P3)\ 
of the relevant half of Definition \3.1\ at xo )■ 

In particular, if u and v are finite for all x G Q and t* < t < T , then they are respectively a viscosity 
subsolution and a viscosity supersolution of (|2.12p on this time interval. 

Proof. When xq G f2, since we can find in O a ((-neighborhood of Xo, the proof follows from [2l|, Propo- 
sition 3.3]. Therefore we shall focus on the case when xo G dil. We give the proof for u. The argument 
for v is entirely parallel, relying on Proposition 14.51 Consider a smooth function (ft such that u — (ft has 
a local maximum at (xo,to). Adding a constant, we can assume m(xo,<o) = <t>(%o,to)- Replacing (ft by 
(ft + \\x — xo\\ 4 + \t — t \ 2 if necessary, we can assume that the local maximum is strict, i.e. that 

(3.3) u(x,t) < <ft(x,t) for < \\(x,t) - (x ,t )\\ < r, 

for some r > 0. By the definition of u, there exist sequences yk, tk = T — NkS 2 such that 

jjk^xo, t k ^to, u Eh (y k ,t k ) -> u(x ,t ). 
Let yk and tk = T — Nks\ satisfying 

(u Ek -<t>){y k ,t k ) > sup [u Ek -(ft){x,t) -e|. 

1 1 (ar,t) — (ac ,*o ) II <T" 
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Notice that since u Ek is defined only at discrete times, the sup is taken only over such times. Evidently, 

(u £k - 4>){Vk,t k ) > {u £k - (f>)(y k> t k ) - 4 
and the right-hand side tends to as k — >• +00. It follows using (|3.3[) that 

(yk, t k ) -> (j/o, *o) and (j/ fc , t k ) -> u(a:o, to), 

as fc — > +00. Setting = u 6k {y k ,t k ) — <j){y k ,t k ), we also have by construction that 
(3.4) £ fc -> and u £k (x,t) < 4>{x,t) + £ k + e\ whenever t = T- n k e\ and ||(x,t) - (x ,t Q )\\ < r. 
Now we use the dynamic programming inequality (|2.25[) at (y k ,t k ), which can be written concisely as 



s {y k ,t k ) < supinf {u £k (y k + Ax,t k + 4) - Az } 
P ,r A * 



with the convention 



Az = p ■ Ax + ^(rAx, Ax) + elf{t k ,y k , u £k (y k ,t k ),p, T) -\\Ax- Ax\\h(y k + Ax). 

Using the definition of £ k , (13. 4p . and the fact that Ax is bounded by a positive power of e, we conclude 
that 

(3.5) 4>(y k ,t k ) + Cfe < sup inf U(y k + Ax, t k + e 2 ) + & + e 3 k - Az} , 

when k is sufficiently large. Dropping £/. from both sides of ()3.5|) , we deduce, by introducing the operator 
S e defined by gj}, that 

(3.6) Mfe.tfc) < S e [y k ,t k ,u £k (y k ,t k ),cf>(;t k +el)}+o(el). 

According to the consistency estimates provided by Proposition ^. 1Q[ we shall introduce four sets (Aj)x<i<4 
respectively defined by 

M := (k G N : d(y k ) < e l k ~ a and M£[<f>(.,t k + $)] > | ||£> 2 ^ fc ,t fe + 4)114"" 
A 2 := |fc G N : 4~ Q "4 < <*(!/*) < and Af£[0(-,t fc +4)] < *fc + 4)H~ 

|j{fceN:% fe ) >4r Q |, 

4:={i£N: d(y fc ) < 4" a - 4 and - e 1 ^ < M^(-,t k + e 2 )] < ^\\D 2 <t>(y k ,t k + 4)||4" 

A, := {fc G N : d(y fc ) < 4"° - 4 and A*™ [<£(•, t fc + 4-)] < -4~"~ K } . 

where p and k are defined in Section l4.1.2l bv (|4.23[) and (|4.24|) and satisfy < re < 1 — a < p < 1. Since 
Ui<i<4v4i = N, at least one set among Ax, A 2 , A 3 and A4 is necessarily unbounded. We shall consider 
these four cases. 

• If A\ is unbounded, up to a subsequence, we may assume that A\ = N. Taking the limit k — > +00, 
we deduce that liminf M^ k [</>(■, t k + e 2 k )] > 0. Since M|* [<f>(-, t k + e 2 k )] -> h(x Q ) - D(f>(x , to) ■ n(x ) by 

k — > + oo 

Lemma l3.3( it follows in the limit k — > 00 that 

(3.7) D(f>(x Q ,t Q ) ■ n(x ) - h{x ) < 0. 

We can notice this result also holds through ()3.6|l . We can apply the second alternative given by (|4.60p 
in Proposition 14. 161 to evaluate the right-hand side of p.6[) . This gives 

<t>(y k ,t k ) ~ 4>{Vk,tk+e 2 k ) < 3el- a Mg[<K;t k +4)] + Ce 2 (l + \u £k (y k ,t k )\) + o(e 2 ), 
where C depends only on and \\D<fi(-,t k + 4)llc 1 (f2ns(t/ fc e 1- "))' Since f° r ^ large enough, 

\\D<P(;t k + 4)\\cl(UnB(y h , s l-")) < SU P ll- D ^(-i i )llci(nnfl(xo,r))> 

\t — tn\<r 
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we can suppose that C depends only on and this sup, which is finite (since <p is smooth) and 

independent of k. By smoothness of <p we have 

-eldt^yuM) +o(4) - C(l + \u*»(y k ,tk)\)4 < 3el- a Mg[<K;tk + 4)]- 
By dividing by e\~ a we obtain 

-el+^dt^h) - C(l + |u e *(lfc.**)l)) + o(el +a ) < 3My*[<f>(;t k +el)}. 

The sequences (u Ek (yk,tk))keK and {dt&iVk, ifc))feeN are respectively bounded by definition of u(xo,to) 
and smoothness of <p. By passing to the limit on k, liminf M^[4>(-,tk + £ 2 )] > 0. By Lemma T3. 31 we 

k— >+oo k 

know that M^[<f)(-,tk + s 2 .)] — > h(xo) — Dcj)(xo,to) ■ n(xo) and (|3.7p is retrieved. 

If A2 is unbounded, up to a subsequence, we may assume that A2 — N. We can apply Proposition ^. 101 



(ii) to evaluate the right-hand side of (|3.6p . This gives 
(t>(Vk,tk) < <t>{yk,tk + sl) ~ £kf( t k,yk,u Ek (y k ,tk),D<j)(yk,tk+£l),D' 2 (l)(yk,tk+£l))+o(el). 
By smoothness of <fi and Lipschitz continuity of / with respect to p and T, we obtain 

Hyk,t k ) - (t>(y k ,tk+et) < -elf(tk,yk,u £k (y k ,tk) : D(t>(yk,tk),D 2 <j)(yk,tk))+o(el). 
It follows in the limit k — > 00 that 

(3.8) d t (i>(x ,to) - f(t , x , u(x ,t ), D<f)(x , t ), D 2 (f>(x , t )) > 0. 

• If A3 is unbounded, up to a subsequence, we may assume that A3 = N. By passing to the limit 

on k, we have that M/ fc *" fe+ k [<f>] tends to zero when e k tends to zero. Since M^ k [4>(-,tk + e 2 ,)] — » 
h{xo) — D(f>(xo, to) -n(xo) by Lemma [3.3[ it follows in the limit k — > 00 that D(j)(xo, to)-n(xo) — h(xo) = 0. 

• If A4 is unbounded, up to a subsequence, we may assume that A4 = N. Hence, taking the limit 
k — > +00, we have 

(3.9) lim sup Mil t^O. ** + £ t)] < 0- 

k— >+oo 

Moreover, by applying the fourth alternative given by (|4. 60[) in Proposition 14.161 to evaluate the right- 
hand side of (13.61). we obtain 



Hvk,t k ) < <f>(y k ,t k + 4) + \4 MV e k M-^ + 4)\ + Ce 2 k (l + \u Sk (y k ,t k )\) + o(4), 
where C depends only on and sup \\D4 ) {',t)\\c 1 (nnB(x r)) * ne same arguments used above 

\t — tc\<r 

for A\. By smoothness of <p we have 

-eldt^ykM) +o(4) - C(l + |« B *(ifo,t*)l)e* < ^M^[0(.,t, +4)]. 
By dividing by e£ we get 

-e2- p (M!/fc.<fc)-C(l + l« e *(w.**)l)) + °(4~ P ) < 

The sequences (u Sk (yk,tk))k& and {dt4>(yk,tk))k^n are respectively bounded by definition of U(xo,io) 
and smoothness of 0. By passing to the limit as k — > +00, we have 

liminfAQW,t fc + e£)] >0. 

k— >+oc 

By comparing this inequality with (|3.9p and using Lemma 13.31 we deduce that 

Dcj)(x ,t ) ■ n(x ) - h(x ) = 0. 

Moreover, we can also apply Lemma [4.91 since e\~ a <C £i~ a ~ K . By the same manipulations as those 
done for the set A2, the inequality (|3.8p holds also true. 
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Thus u is a viscosity subsolution at (xo,to). 



We turn now to case (iii) i.e. the case to = T. If xo G Ct, the analysis led in [2lj, Proposition 3.3] gives 
the result. It remains to study u on the boundary. We want to show that u(-,T) = g is also satisfied 
on dCl. By the definition of u given by (|2.27[) and considering a particular sequence (ek,Xk,tk — T)k<£N 
which converges to (0, xo, T), it is clear that u(-, T) > g on dfl (using the continuity of g and the fact that 
each u £ has final value g). If this sequence realizes the sup, we have in fact the equality. The preceding 
argument can still be used provided tk < T for all sufficiently large fc. Thus, considering the different 
possibilities according to tk < T or tk — T and also on Xk G fi or Xk G dQ, we know that for any smooth 
4> such that u — 4> has a local maximum at (xq, T), 

(3.10) either u(xq,T) = g(xo) or else 

max (dt<p(x ,T) - f(t , x , u(x , T), D(/)(x , T), D 2 4>(x , T)), h(x ) - Dcf)(x , T) ■ n(x )) > 0. 

Moreover this statement applies not only at the given point Xq, but also at any point nearby. 

Now consider the functions 

u 4\ -i +\ 11^ ~ x oll 2 T ~ t , V AI \ 
ip[x, t) — u(x, t) h Ka[x) 

T) [I 

and 

(3.11) 4>{x, t) = I|X ~ X ° 112 + — - - Kd(x), 

where the parameters 77, /1 are small and positive and K — ||/i||z,oo + 1. Suppose u is uniformly bounded 
on the closed half-ball {||(x, t) — (xo,T)\\ < r,t < T} and let tp attain its maximum on this half-ball at 
(x v .n,t v ,fj,)- We assume r is small enough such that d is C 2 on this half-ball so that <j> can be taken as a 
test function. We clearly have 

(3.12) u(x v ^,t, hfl ) + Kd{x, q ^) > ip(x v ^, t V!ll ) > ip(x , T) = u(x , T). 

By plugging the expression of ^(iKjj,^, in the right-hand side of inequality (|3.I2I) . we obtain 

(3.13) 0< W x ™- x °W 2 + T - t *>» <u(x rh ^t r ^)~u(x Ql T) + Kd(x^). 

V M 

Since u is bounded on the half-ball and belongs to the half ball for all 77 and /i, the right-hand side 
of (|3.13[) is bounded independently of 77, /j, which yields 

(3.14) {x w ,t w ) ->• (x ,T) as 77, ft->0. 

By using the upper semicontinuity of u and taking the limit (13.1411 in (|3.12j) , we get 

(3.15) u(x Vjfl ,t V)IJ ,) ->• u(xq,T) as?7,/^^0. 
By combining (|3.14D and (|3.15[) . we finally obtain by (|3 . 13|) that 

(3.16) W»-xo\\ 2 + T-t Iht ^ Q a 

77 /j 

If t v ^ < T and x n ^ G Cl then part (i) of Proposition 13.41 applied to 4> defined by (|3.11[) assures us that 

(3.17) x Vjtl , u(x v ^, t m „), 2 ^ - X ° - KDd(x^), —I — KD 2 d(x v ^)) > 0. 

[A V V 

Since / is continuous, for any 77 > there exists /i > such that (|3.17[) cannot happen. Restricting our 
attention to such choices of 77 and /j, it remains to examine two situations: on the one hand t Vjtl < T and 
Xr h n G dQ and on the other hand ^ = T. Arguing by contradiction, let us assume that t^^ < T and 
x ri,n G dQ. By the Taylor expansion on the distance function close to xq, we have 

d(x) = d(x ) + Dd(x ) ■ (x — x ) + 0(\\x — x \\ 2 )- 
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By using that xq and x v ^ are on the boundary <90, d(xo) = d(x r) ^) = and Dd(xo) = — n(xo), this 
relation reduces to 

(3.18) n(x ) ■ {x Vtli - xq) = 0(\\x v ^ - x \\ 2 ). 

By combining (|3.16p and f|3 . 18[) . we compute 

2 

D4>(x v> n,tq tll ) ■ n(x ) = -(x VtM - x Q ) ■ n(x ) - KDd(x Vlll ) ■ n(x ) 

= O OjElhiiZJ^ \ + Kn{x n ,n) ■ n(x ) -> K , as rj, fi -> 0. 

By smoothness of 4> and continuity of n on SO, we deduce that Dcp^x^^, T)-n{x ri>ll ) — > ||/i||l°°+1 > M 2 ^,/*) 
which denies the second alternative proposed at (|3.10[) . As a result, the only remaining possibility for 
(|3.10p is u(x v>IJl ,T) — g(x r] ^). By continuity of g, it follows in the limit 77, fi — > that u(xq,T) = g(xo), 
as asserted. □ 



Remark 3.5. In the proof of the convergence at the final-time in Theorem \3.4\ we needed in a essential 
way that the domain was assumed to be at least C 2 . More precisely, in this case, since the distance 
function d is C 2 in a neighborhood of the boundary, it allows us to take tf> given by (|3 . 1 1 1) as a test 
function. 



3.3. The elliptic case. We turn now to the stationary setting discussed in Section 1231 As in the time- 
dependent setting, our convergence result depends on the fundamental consistency result Proposition 
14.171 So we require that the parameters a, /?, 7 satisfy (|2.13p - (|2.15p . and that f(x,z,p,T) satisfy not 
only the monotonicity condition (|1.4p but also the Lipschitz continuity and growth conditions (jl.9p - 
(jl.lOp . To prove that U e is well defined, we require that the interest rate A be large enough, condition 
(jl.8p . and that h be uniformly bounded. Finally, concerning the parameters m and M and the function 
tjj associated to the termination of the game, we assume that ip € C 2 (0) fulfills §^ = \\h\\ L -> + 1 on 50, 
m = M — 1 — 2|| , 0|| Loo ^, x = m + 1 1 ip II £00 (jf) + ip and M is sufficiently large. These hypotheses ensure 
us the availability of the dynamic programming inequalities stated in Proposition [ 



Proposition 3.6. Suppose f , g, A and a, f3, 7, m, M , -0 satisfy the hypotheses just listed (from which 
it follows that v and u are bounded away from ±M for all i£f!j. Then u is a viscosity subsolution and 
v is a viscosity supersolution of (I2.54p in Q. More specifically: 



if Xq G fl then each of u and v satisfies part (El ) of relevant half of Definition \3.2\ at xq, and 



if xq € dfl then each of u and v satisfies part (E2) of relevant half of DeHnition \3.S\ at xq. 



Proof. When xq € f2, the proof is similar to that of Theorem 13.41 Therefore we shall focus on the case 
when xq £ dfl. We give the proof for u, the arguments for v being similar and even easier due to fewer 
cases to distinguish. Consider a smooth function <fi such that u— <f> has local maximum on 17 at xq € dQ. 
We may assume that (D(f>(xo),n(xQ)) > h(xo) since otherwise the assertion is trivial. Adjusting 4> if 
necessary, we can assume that u(xq) — 4>{xq) and that the local maximum is strict, i.e. 

(3.19) u(x)<(f>(x) for x £ On {0 < ||x- z || < r}, 

for some r > 0. By the definition of It given by (|2.65p . there exist sequences > and jjk G O such that 

i/k^-XQ, u Ek (y k ) u(x ). 
We may choose y k G O such that (u ek — 4>){yk) > sup (u £k — 4>){ x ) — £ k- Evidently 

fTn{||2;-2;o||<r} 

(u £k ~4>){y k ) > (u £k - 4>){m) - 4 

and the right-hand side tends to as k —> 00. It follows using (|3.19p that 

y k -t x Q and u ek (y k ) -t u(x ), 
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as k — > oo. Setting = (u £k — (j))(yk), we also have by construction that 

(3.20) & and u £h (x) < <f>(x) + £fc - eg whenever ieS] with ||x - x || < r. 

We now use the dynamic programming inequality (j2.62[) at y k , which can be written concisely as 

u £k (y k ) < supinf {e- X£k u £k (y k + Ax) - S k \ , 

p p Ax I J 

with the convention 

S k =p - Ax + ~{TAx,Ax) + e 2 k f(x,u £k {x),p,T) - \\Ax - Ax\\h(x + Ax). 

By the rule f|2 . 18|) of the game, for every move Ax decided by Mark, the point y k + Ax belongs to Q. 
Combining this observation with (|3.20p and the definition of we conclude that 

4>{Vk) + & < supinf \e~ x ^ [4>{y k + Ax) + £ fc - eg] - 6 k \ . 

p p Ax I ) 

We may replace e _A£fc by 1 — \e\ and e~ X£k £ k by while making an error which is only o(e 2 ) using the 
fact that £fc — > 0. Dropping £ k from both sides, we conclude that 



4>(Vk) < supinf (e- XEk (j)(y k + Ax) - 5 k ) + o(e 2 k ). 

p P Ax \ / 



We can evaluate the right-hand side using Proposition 14.171 case (ii) for k large enough. We introduce p 
and k defined in Section T4.1.2I by (|4.23|) and (|4.24[) and satisfying < K < 1 — a < p < 1. If we may 
assume, up to a subsequence, that for all k large enough, on the one hand d(y k ) > £ k ~ a or on the one 
hand e\~ a - e p k < d(y k ) < e\~ a and Mg[(/>] < ^\\D 2 <t){y k )\\e k ~ a , we can apply Proposition \ZJT\ case p) 
to evaluate the right-hand side 

< u £k (y k ),D(f>(y k ), D 2 cp{y k )) - e 2 \u £k (y k ) + o(e 2 ). 

By passing to the limit k — > +oo, we get the required inequality in the viscosity sense. Otherwise, recall 
that (D(j)(xo),n(xo)) > h(xo). By Lemma l3.3i we have 

(3.21) Mf k W\ h ( x o) ~ (D<Kxo), n(x )) < 0, 
and the condition 

(3.22) M£[cf>]<-e 1 k - a - K 

is satisfied for all k sufficiently large. Therefore, up to a subsequence, it remains to consider a sequence 
(y k , £fc)/ceN satisfying both d(y k ) < s k ~ a — e k and (|3.22[) . The last part of Proposition 14. 1 81 can be applied 
and we get by (|4.7ip that there exists a constant C depending only on M, \\D(f>\\ c i,££, nB , yh E i-«s and 
||/i||loo such that 

< \e p k MH [0] + Ce 2 - \e 2 <jy{y k ) + o(e 2 ), 
recalling that (e 1 ^" — d(y k )) > e k and MJ** [<p] < 0. By dividing by e kl it follows that 

-e 2 k - p (C-my k ))+o(e 2 k - p ) < jMj£[# 

The sequence {4>{yk))ken is bounded by smoothness of <f>. Since \\D<l>\\c 1 ((i)nB(v h ,e 1 - a ) - \\D < f'\\c 1 (n)nB(xo,r) 
holds for k large enough, we can assume that C is independent of k depending only on ll-D^llc^njnB^o r)> 
M and ||/i||j,oo. Taking the limit as k — > +oo, we deduce that 

liminf M Vk U\ > 0, 

which is a contradiction with (|3.2ip . Thus u is a viscosity subsolution at xq. □ 
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4. Consistency 



A numerical scheme is said to be consistent if every smooth solution of the PDE satisfies it modulo an 
error that tends to zero with the step size. It is the idea of the argument used in [2l|. In our case, we must 
understand how the Neumann condition interferes with the estimates proposed in [U Section 4]. The 
essence of our formal argument in Section 12.2.11 was that the Neumann condition term is predominant 
compared to the PDE term at the boundary and produces a degeneracy in the consistency estimate. The 
present section clarifies the connection between our formal argument and the consistency of the game, 
by discussing consistency in more conventional terms. The main point is presented in Propositions 14.51 
and 14.101 In order to explain very precisely how the consistency estimate of (ill Section 4] degenerates, 
we establish the consistency of our game as a numerical scheme by focusing on different cases according 
to the values of the quantities m*[fl and M^[c/>] defined by (|3.1|l - (|3.2|l and the distance d(x) to the 
boundary d£l. 



4.1. The parabolic case. Consider the game discussed in Section l2.2l for solving —u* + f(t. x. u. Du. D 2 u) = 
in £1 with final-time data u{x,T) = g(x) for x G £1 and boundary condition §£(x,t) = h(x) for 
x 6 dtt,t € (0,T). The dynamic programming principles (|2.25l) - (|2.26p can be written as 



where S e [x, t, z, <fi] is defined for any x G fi, z € M. and t <T and any continuous function <j>: Q, — > M by 
(4.1) S £ [x, t, z, (j)] = maxmin [(f> (x + Ax) 



subject to the usual constraints ||p|| < e °, \\T\\ < e 7 , || Ax|| < e 1 a and Ax = projjy(x + Ax) — x. The 
operator S e clearly satisfies the three following properties: 

• For all 4> £ C(U), S [x, t, z, 0] = (j>{x). 

• S e is monotone, i.e. if tf>i < 02, then S e [x,t,z,<f>i] < S £ [x,t,z,(f>2\- 

• For all (j) G C(H) and c6l, 

(4.2) S e [x, t,z,c + <f>] = c + S e [x, t, z, <f)] . 

Fixing x, t, z and a smooth function <f>, a Taylor expansion shows that for any || Ax\\ < e 1 ~ a , 



u £ (x,t) < S e [x,t 7 u £ (x,t),u £ (-,t + e 2 )] , 
v £ (x,t) > S e [x,t,v £ {x,t),v £ (- 1 t + e 2 )] , 



p,r Ax 



(• 




') 




Since a < 1/3 by hypothesis, e 



.3 — 3q 



o(e 2 ). By rearranging the terms, we compute 



4> (x + Ax) - 



p- Ax + - (TAx,Ax) +e 2 f (t,x,z,p,T) - \\Ax- Ax\\h(x + Ax)j 
<f>(x) + (D(f>(x) - p) ■ Ax + \\Ax - Ax|| [h(x + Ax) - D(j)(x) ■ n(x + Ax)} 

+ X - {D 2 (b(x)Ax, Ax) - i (TAx, Ax) - e 2 f(t, x, z,p, V) + o(e 2 ), 
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Aa; — Ax 

since the outward normal can be expressed by n(x + Ax) = — — — - if x + Aa; 4 Q and the move 

II Ax — Ax\\ 

Ax can be decomposed as Ax = Ax + (Ax — Ax). Thus, we shall examine 



(4.3) S £ [x, t, z, 4>] — 4>{x) — maxmin 

p,r Ax 



{Dcf>(x) -p)-Ax+- (D 2 (j>(x)Ax, Ax) 



+ ||A1- Ax\\{h(x + Ax) - D(f>(x) ■ n(x + Ax)} - - (TAx,Ax) - e 2 f(t, x, z,p,T) 



+ o(s 2 ). 



4.1.1. Preliminary geometric lemmas. This subsection is devoted to some geometric properties of the 
game that will be useful to show consistency in Section 14.1.21 We start by some estimates, involving the 
geometric conditions on the domain, about the moves Ax decided by Mark. 

Lemma 4.1. Suppose that ft is a C 2 -domain satisfying the uniform exterior ball condition for a certain 
r > 0. Then, for all < e < r 1 - a and for all Ax constrained by (|2.17j) . determining Ax by (|2.22ll . we 
have 

(4.4) ||Ax-Ax|| ^e 1 -"-^) and || Ax|| < 2e 1 - a - d{x). 

Proof. Let us prove the first inequality, the second following immediately by the triangle inequality. If 
the point x = x + Ax belongs to Q, Ax = Ax and the result is obvious. Supposing now x does not belong 
to fl, the set S = [x, x] n <90 is not empty and we can consider a point xq £ S. By the rule of the game, 
we have ||x — x|| = ||Ax|| < e 1 "". Since xa £ dfl by construction, it is clear that ||x — xo\\ > d(x). We 
deduce that 

\\xg - x\\ = \\x - x\\ - \\x d - x\\ < e x - a - d(x). 

By the uniform exterior ball condition, the orthogonal projection on is well-defined on fl(e 1 ^ a ) C ^i(r). 
By property of the orthogonal projection and since x ^ fi, we can write 

|| Ax — Ax|| = mfjly — x\\ = inf \\y — x\\ < \\xg — x\\, 
yen yedn 

which gives directly the first estimate in (j4.4[) . □ 

The following lemma uses the crucial geometric fact that £1 satisfies the interior ball condition intro- 
duced in Definition 11.21 for which there is no neck pitching for e sufficiently small. 

Lemma 4.2. Let a > I — a and B > 0. Suppose that fl is a domain with C 2 -boundary dil and satisfies the 
uniform interior ball condition. Then, for all possible moves \\Ax\\ < e 1 "" such that \\ Ax + e 1 " Q! n(x)|| < 
Be a the intermediate point x belongs to Q for all e sufficiently small. Moreover, for all possible moves 
||Ax|| < such that \\Ax — e 1_Q n(x)|| < Be a and Ax determined by (|2.22[) , we have 

(4.5) || Ax - Ax|| > e 1 -" - d{x) - Be a + 0(e 2 ~ 2Q ). 

Furthermore, if in addition we assume d(x) > £ 1_Q — e n with 1 — a < rj < a , the intermediate point x is 
outside for all e sufficiently small. 

Proof. For the first assertion, since D, satisfies the uniform interior ball condition (there is no neck pitching 
for e sufficiently small), we observe that the set <9f2 n B(x, 2e 1 ~ a ) is below a paraboloid Pi of opening A 
and above a paraboloid P-2 of opening —A touching dil at x. By the Taylor expansion, if T s dVl denotes 
the tangent space to dil at x, we get that for all y £ d£l (1 B(x, 2e 1_a ), 

\(y-x)-n(x)\=d(y,T s dQ) < ^A{2e 1 - a ) 2 1 

Since (x + Ax — x) • n(x) < — — d(x) + Be a , we deduce that for all e sufficiently small, 

(x + Ax — x) • n(x) < inf (y — x) ■ n(x), 

yednnBix^e 1 —) 



which yields that x + Ax belongs to Q. 
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For the second claim, we denote by (Kj(x))i<j<jv-i the principal curvatures at x on dCl and by 
(ei, • • • , ejv) an orthonormal frame centered in x with first vector e\ — n(x). Since CI is a C 2 -domain, 
(ea, • • • , ejv) form a basis of the tangent space T x dCl. We compute 

e 1 ^ - Be a < Ax ■ n(x) = (Ax - e 1 ""^)) • n(x) + 

Thus x is contained in the half-space Hi determined by (y-x)-ei > e 1 ~ a — d(x) — Be a and d(x, T^dCl) > 
e 1_Q — d(x) — Be" . Moreover, we deduce from (|4.4[) and the triangle inequality that for each move Ax 
we have x + Ax G B(x, 2e 1 ~ Q ). Assume xi — p(x2, ■ • • , xn) is a local C 2 -parametrization of dCl around 
x. By a Taylor argument and by continuity of the principal curvatures on dCl, it follows that, for e > 
small enough, 

(4.6) d(x + Ax, T s dCl) < i(7 1 (2e 1 - Q ) 2 = 2C ie 2 - 2Q , 

where Ci := 2 max { | (x) | : I < i < N — 1}. By the triangle inequality, we deduce that 
\\x + Ax — x|| > ||proj T _ ao (x + Ax) — x|| — ||x + Ax — proj T _ aa (x + Ax)|| 

> d(x, T^dCl) - d(x + Ax, T^dCl) 

> e 1 "" - d(x) - Be a - 2C l£ 2 - 2a . 

In particular, if d(x) > e 1_Q — e'' with 1 — a < rj < a the right-hand side is strictly positive for e 
sufficiently small and x ^ Cl. □ 

The next lemmas gather some estimates which will be useful to establish our consistency estimates. 

Lemma 4.3. Under the hypothesis of Lemma \J h .2\ for all moves Ax constrained by (|2.17|l . determining 
Ax by (|2~22l . we have 



(4.7) - -(s'- a - d(x)) < -- I 1 - 1 (Ax) • n(x) + ||Ax - Ax|| < -(e'- a - d(x)). 



Proof. The left-hand side of (|4.7p can be written in the form 

l(Ax)-n(x) ||Ax-Ax|| 



2 e 1 "" e l ~ a - d(x) 



which directly gives the desired estimates by using (|2.17p and the first inequality given by (|4.4[) . □ 

Lemma 4.4. Let A G M JV (R), k G C b (dCl) extended by some function k : Cl — » M, and x G f2. Suppose 
in addition that 

(4.8) (3e 1 -"-rf(x))||A|| < inf fc(x + Ax), 

A3 

im£/i Ax constrained by (|2.17p and Ax determined by (|2.22p . TTien 

(4.9) min{(ylAx, Ax) - (AAx, Ax) + ||Ax - Ax\\k(x + Ax)} = 0, 

Ax 

where Ax is constrained by (|2.17p and determines Ax by (|2.22p . 

Proof. If x = x + Ax G fi, the function is equal to zero. We now consider the moves for which x ^ Cl. 
Then 

(4.10) (AAx, Ax) - (AAx, Ax) = (A(Ax - Ax), Ax - Ax) + 2 (AAx, Ax - Ax) . 
By the Cauchy-Schwarz inequality, we obtain 

(4.11) | (AAx, Ax) - (AAx, Ax) \ < \\ A\\ \\ Ax - Ax\\ (||Ax- Ax|| + 2||Ax||). 
By using (|4.4p and ||Ax|| < e 1_Q , we get 

(4.12) | (AAx, Ax) - (AAx, Ax) \ < \\A\\\\Ax ~ Ax\\ (Se 1 "" - d(x)) . 
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Thus 



(AAx, Ax)-(AAx, Ax) + \\Ax-Ax\\k(x+Ax) > \\Ax-Ax\\ { inf k(x + Ax) - ||A||(3£ 1_Q - d(x)) 

Ax 

The right-hand side of this last inequality is strictly positive by the assumption (|4.8j) . □ 



4.1.2. Consistency estimates. In this subsection we state our consistency estimates. They explain pre- 
cisely the conditions under which the usual estimate proposed in |2ll | holds for x near the boundary and 
<f» £ C 2 (£l). If it does not hold, there is a degeneration of the estimates respecting the final discussion of 
formal derivation of the PDE at Section [2. 2. II For fixed x £ Slfe 1- "), these estimates take into account 
the size and the sign of the boundary condition in the small ball B(x, e 1_Q ) and the distance d(x) to the 
boundary. 

In the heuristic derivation presented in Section 12.2.11 we assumed that Ax i— > h{x + Ax) — D<p(x) ■ 
n(x + Ax), with Aa; determined by (|2.22j) . was locally constant in a ^-neighborhood of the boundary 
near x. In the general case, this hypothesis must be relaxed. To do this, we observe that, for all Ax 
constrained by (|2.17[) satisfying x + Ax ^ fl and determining Ax by (12.22)) , 

m%[<t>] < h(x + Ax) - D(j){x) ■ n(x + Ax) < [(f)] , 

where m*[cf)] and M*[4>] are defined by (13.1I) - (|3.2I) . Therefore we are going to specify some strategies for 
Helen which are associated to the two extreme situations mff [(f)] and [(f)] by following the optimal choices 
(|2.48l) and (|2.52| obtained in the formal derivation at Section f2. 2. II More precisely, for all x £ U(e 1 ~ a ), 
we define the strategies p™ t (x), p%^ t (x) and r opt (x) in an orthonormal basis 23 = (ei = n(x), e%, ■ ■ ■ , ejv) 
respectively by 



(4.13) 

(4.14) 

and 

(4.15) 



p™ t (x) = D(f>(x) + 



p% t (x) = D(f>{x) 



1 



1 - 



d(x) 
d(x) 



rrfM- 

mm - 



4 

-1-Q 



1 - 



1 - 



d\x) 

£ 2-2a 

d 2 {x) 

c 2-2a 



(D'(f>(x))u 
{D 2 (f)(x)) u 



n(x), 
n(x), 



opt 



(x) - D 2 (f>{x) + 



; -i 



d 2 (x) 



(D 2 (f>(x))\ 



E 



U) 



where En denotes the unit-matrix (1,1) in the basis 23. These strategies depend on the local behavior 
of (f> (size and amplitude) around the boundary and on the geometry of the boundary itself. 

Since there is a degeneration of the usual estimates, there is no hope for one simple estimate. We are 
going to separate the study in two steps: Proposition 14.51 provides the estimates for the lower bound and 
Proposition 14. 101 deals with the upper bound. Moreover, Section [4~2l is devoted to the technical proof of 
the upper bound distinguishing several cases according to the size of Mff: [(f)] and d{x). 



Proposition 4.5. Let f satisfy (|1.4p and (|1.6p - (|1.7p and assume a, f3, 7 satisfy (|2 . 13[) (|2 . 1 5[) . Let p™ pt 
and r pt be respectively defined in the orthonormal basis (ei = n(x), e2, • • • , e^) by (|4. 1 3[) and (|4.15p . For 
any x, t, z and any smooth function (f> defined nearx, S e [x,t, z,(p] being defined by (|4. 1|) . we distinguish 
two cases: 

i. Big bonus: if d(x) > e 1 "" or m x e [(f)] > i^e 1 "" - d(x))\\D 2 (f>(x)\\, then 

S e [x, t, z, 4>] - (f)(x) > -e 2 f(t, x, z, D(f>(x), D 2 (f){x)). 

ii. Penalty or small bonus: if d(x) < e 1- " and m^[(p] < ^(3e 1_Q! — d(a;))||£' 2 (/)(x)||, then 

S s [x, t, z, (f)] - (f>{x) > ^(e 1 -" - d(x)) (sm x e [(f>] - A\\D 2 ^{x)^ < 
where s = — 1 if mfff[(f>] > and s = 3 if m^[(f)] < 0. 



») - e 2 /(i, x, z,p™ t (x), T opt (x)), 
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Proof. If d(x) > e 1 a , the usual estimate [2ll . Lemma 4.1] holds. We now focus on the case d(x) < e 1 a . 
By the definition of m x [(f)] given by ()3.1j) and the positivity of ||Ax — Ax\\, for all ||Ax|| < e 1_Q , we have 

(4.16) ||Ai - Ax\\ {h(x + Ax) - D<j>(x) ■ n(x + Ax)} > \\Ax - Ax\\m*[<j>]. 

Therefore it is sufficient to find a lower bound for 



max mm 

p.r Ai 



(Dcj)(x) -p) ■ Ax + m x [(f>][\Ax - Ax[\ + )^{D 2 <j>(x) Ax , Ax) ~ ^(TAx,Ax) - e 2 f(t,x, z,p,T) 



where p, T and Ax are constrained by (|2.16[) - (12.17p and Aa; determined by (|2.22[) . In other words, by 
taking advantage of the monotonicity of the operator S e with (|4.16p , we shall look for a lower bound for 
an approximated operator bounding S e from below and very close to it when e — > 0. 

Then, we also observe that for every choice p and T, 



S e [x,t,z,4>] - <f>(x) > ~e 2 f(t,x,z,p,T) 
+ min 



(Dcj)(x) - p) ■ Ax + i (D 2 cj)(x)Ax 7 Ax) - ^ (rAx, Ax) + \\Ax ~ Ax\\ 



We now distinguish two particular strategies for Helen. For part (i) we consider the particular choice 
p = D(j){x), T = D 2 (f)(x) and obtain 

S e [x, t, z, 0] - <f)(x) > -e 2 f (t, x, z, Dcj)(x),D 2 <j)(x)) 



mm 

Ax 



h(D 2 cf)(x)Ax, As) - i (D 2 cj)(x)Ax, Ax) + \\Ax- Ax\\m* [</>] 



> -e 2 f(t,x,z,D<t>{x),D 2 cf>{x)) , 



by applying Lemma l4~4l with A = ^D 2 (f>(x). For part (ii) we consider the choice p = p™ t (x), T — r opt (a;) 
and find 



x) > -e 2 f(t,x,z,p™(x),T opt (x))+l x [4,], 



with l x [4>] defined by 
(4.17) 



l x [(b] = min 
Ai 



(Dcf)(x) — p™ t ) • (Ai) + \ (D 2 (f>(x)Ax,Ax) - \ (T opt (x)Ax, Ax) + \\Ax - Ax\\m x [4>] 



It now remains to give a lower bound for l x [(f)] . By plugging the expression f|4. 13[) of t (x) in (|4.17l) , we 
have 



mm 

Ax 



-1-a 



d(x) 



2e 



l-a 



(Ax) 



\Ax - Ax\\ m x [(j>} 



+- (D 2 (j,(x)Ax,Ax) - - (T opt (x)Ax,Ax) + 



d 2 {x) 



(i? 2 </ ) (x)) 11 (Ax) 1 



It is clear that l x [(f)] > If [(f)] + l x [(f)] with If [(f)] and l x [(f)] respectively defined by 



(4.18) 
and 

(4.19) lj[fl:=~mia 

Z Ax 



mm 

Ai 



d(x) 



2s 1 



(D 2 (f>(x)Ax,Ax) 



(Ax)! + || - Aa:|| ) m x [(f)] 
d 2 (x) 



opt 



(x)Ax,Ax) 



1 



(D 2 ^(x))n(Ax), 



By using Lemmas 14.61 and 14.81 stated below, giving lower bounds respectively for If [(f)] and l x [(f>], one 
obtains 

l x [(p] > ^e 1 "" - d(x))m x M - 2\\D 2 (f>(x)\\e 1 - a (e 1 ^ - d(x)) = ±(e l ~ a - d(x)) (sm x [(f>] - 4||# 2 #r)|| £ 1 - Q ) , 



which gives the desired estimate. 



□ 
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The three following lemmas provide the required estimates for lf[<f>] and Zf [4>\- 

Lemma 4.6. For any x £ 0(£ 1_a ) and any function (ft defined at x, If [(ft) being defined by (14.1811 , we 
have 

^ie 1 -" - d(x))mZ[4>] <IM <0, 
with s = — 1 ifm^[(ft] is positive and s — 3 if m^ty] is nonpositive. 

Proof. By considering Ax = 0, If [(ft] is negative. To find a lower bound on If [(ft], if m"[i^] is negative, we 
may write 



2s 



J^(Ax) 1 + ||Ax-Ax| 



d(x))m x e [(ft], 



the last inequality being provided by the right-hand side inequality given in Lemma [4 . 31 since by hypothesis 
m x [(ft] is negative. If [<f\ is nonnegative, the result follows from applying the left-hand side inequality 

□ 



given in Lemma 14.31 

Lemma 4.7. Let x E il(e 1 ^ a ) and (ft 6 C 2 (il). For all Ax constrained by (|2.17p . we have 



(4.20) 

and 

(4.21) 



^(D 2 (ft(x)Ax, Ax) - ^(D 2 (ft(x)Ax, Ax) 



1 



< -\\D 2 q>(x)\\ (Se^-dix)) ||Ax-Ax||, 



^(D 2 (ft(x)Ax,Ax) - ^(T opt (x)Ax,Ax) 



< -\\D 2 <b(x)\\ [e 



l-a 



d(x)) (le 1 -* - d(x)) 



where T opt (x) is the optimal choice defined by (|4.15p in an orthonormal basis 23 = (ei = n(x), ■ ■ ■ , ejv)- 

Proof. The first inequality is an immediate consequence of (|4.11|) . For the second inequality, all the 
coordinates {(D 2 (ft(x) — T opt (x))ei,ej) in the basis 23 are equal to zero, except for i = j = 1. By using 
the vector decomposition given by (|4.10|) . we have 

^{D 2 (ft(x)Ax,Ax) - i(r opt (x)Ax,Ax) = ^{D 2 (ft(x) - r opt (x))n|(Ax)i| 2 

+ -\\Ax - Ax\\ 2 ((D 2 (j)(x)n(x + Ax),n(x + Ax)) - ||Ax - Ax\\{D 2 (ft{x)n{x + Ax), Ax). 



Since {D 2 cft{x) - r opt (x))n = i (l - ^l)(D 2 ^))ii by one obtains 



^(D 2 (ft{x)Ax,Ax) ~ ^(T opt (x)Ax,Ax) 



< \\DU(x) 



1 



d 2 {x) 

c 2-2a 



|(A£)i| 2 + -II Ax - Ax|| 2 + ||Ax - Ax||||Ax| 



The estimate ([4T2T]) now follows from (f2~T7| and ([O 



□ 



Lemma 4.8. For any x € ^(e 1 a ) and any function (ft defined at x, l x [(ft) being defined by (14. 19[) . we 
have 



-2\\D 2 d>(x)\\s 1 - a (e 1 - 



■d(x))<q[<t>]<o. 



Proof. By considering Ax = 0, Z2 is negative. We seek now to find a lower bound on fe- By combining the 
triangle inequality with Lemma l4.7[ the explicit expression of r op t(x) given by ([4~T5| and ||Ax|| < e 1_Q , 



A GAME INTERPRETATION OF THE NEUMANN PROBLEM FOR FULLY NONLINEAR EQUATIONS 



37 



we deduce that 
1 



2 



D 2 cf>(x)Ax, Ax)- (T opt (x)Ax,Ax) + \ (e 1 ^ - ) (£> 2 # > ) , , .( A.D , 



< ill^^IKe 1 -" - d(x)) {le 1 -* - d{x)) + \\\D^(x)\\ (e 1 -" - ^) e 1 ^ 
(4.22) <2\\D 2 cf,(x)\\(e 1 - a -d(x))e 1 ~ a , 

which is precisely the proposed estimate. □ 



We shall now provide the consistency estimates about the upper bound of (|4.3[) . Before stating our 
main estimate in Proposition 14. 101 we can give a simple case for which the usual estimate holds. 

Lemma 4.9. Let f satisfy (|1.4p and (|1.6|) - (|1.7|l and assume a, f3, 7 satisfy (|2.13|) - (|2.15|) . For any x, 

t, z and any smooth function <fi defined near x, S e [x,t, z, <f>] being defined by (|4.1j) . if d(x) < e l ~ a and 
M*[4>] < -±\\D 2 (j){x)\\ (Se 1 - - d(x)), then we have 

S e [x, t, z, (/>} - <p(x) < ~e 2 f(t, x, z, Dcf>(x), D 2 <j)(x)) + o(e 2 ). 

Moreover, the implicit constant in the error term is uniform as x, t and z range over a compact subset 
ofQ, xlxl. 

In the rest of the section, we now accurately focus on the case d(x) < e 1_Q . The goal is to obtain 
precise estimates on (|4.3[) in the following three cases: MJ 5 [0] very negative, M*[<fr] very positive and 
Mj?[(j>] close to zero, the bounds between the cases depending on some powers of e. We have formally 
shown in Section 12.2.11 that the first case is favorable to Mark since Helen can undergo a big penalty if 
Mark chooses to cross the boundary. On the contrary, the second case is preferable to Helen because she 
can receive a big coupon if the boundary is crossed. In the last case, the boundary is transparent (think 
of Mg [(j)] = 0) and the penalization due to the boundary is to be considered only through second order 
terms. In order to establish the precise estimates, we successively introduce two additional parameters 
p, k > such that 

(4.23) 1 - a < p < min ( 1 - 7 ^ ~ ^ , 2 - 2a - 7 



and 

(4.24) j + p-(l-a)<K<l-a. 

These coefficients are well-defined by virtue of (|2.13[) and (|2.15|) . 

Proposition 4.10. Let f satisfy (11.41) and (|1.6|| - (|1.7|) and assume a, (3, 7, p, n satisfy (|2.13p - (|2.15p and 
(|4.23l) - (|4.24l) . Let p^p t and T opt be respectively defined in the orthonormal basis (e% = n(x), e%,- ■■ , cn) 
by (|4.14p and (|4.15p . For any x, t, z and any smooth function (j) defined near x, S e [x, t, z, </>] being defined 
by (|4.ip . we distinguish four cases: 

i. Big bonus: If d(x) < e 1 "" and M*[<j)] > ^\\D 2 (f)(x)\\e 1 - a , then 

S s [x,t,z,4>]-<i>(x) < ^e x - a -d{x))M x e [4>]-e 2 f{t : x 1 z,pl I vt {x),Y opt (x))+o{e 2 ). 

ii. Far from the boundary with a small bonus: if e 1 ~ a —E p < d(x) < e 1_Q and [(f)] < ^\\D 2 (f>(x)\\e 1 ~' 
or if d(x) > e x ~ a , then 

S e [x,t,z,4>] - cf)(x) < -s 2 f(t,x,z,D(f>(x),D 2 <p(x)) +o(e 2 ). 

Hi. Close to the boundary with a small bonus /penalty: if d(x) < e 1 ^" — e p and — e 1 ~ Q_K < M^[<f>] < 



|||U a 0(x)||e 1 - tt , then 



S s [x,t,z,4>] - <f>(x) < -£ jS f(t,x,z,D(t>(x),D' ! 4>(x) + CiI) + o(e 2 ), 

with Cl = f\\D 2 0(x)\\ (l-^iy 
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iv. Close to the boundary with a big penalty: if d(x) < e 1 a — e p and M £[</)] < —e 1 a K , then 

(4.25) S e [x,t,z,cj>]-^(x)< -(e 1 - -d{x))M*[<j)]-e 2 min f(t, x, z,p, T opt {x)) + o(e 2 ), 

4 " peso» t (x),r) 



with r defined by r = - p^l) \M^[(j>]\. 



Moreover, the implicit constants in the error term is uniform as x, t and z range over a compact subset 
ofU xlxl. 

Before proving these estimates, it is worth drawing a parallel with the formal derivation done at 
Section 12.2.11 The lower bound proposed by Proposition 14.51 case (i) corresponds to the formal analysis 
when m > 0. The upper bound proposed by Proposition ^. lOl case (iv) is associated to the formal analysis 
when m < 0. Furthermore, we can observe in the proof that the factor 1/4 in (|4.25p could be replaced 
by any number in [1/4, 1/2), the bound 1/2 corresponding to the heuristic derivation given by (|2.53[) . 

4.2. Proof of Lemma 14. 91 and Proposition 14.101 For sake of notational simplicity, we write A m i n (^l) 
for the smallest eigenvalue of the symmetric matrix A and we omit the x-dependence of p^ t (x) and 
r pt(x). Moreover, by the definition of M*[4>] given by (|3.2p and the positivity of ||Ax — Ax||, for all 
||Ax|| < we have 

(4.26) || As - Ax|| {h(x + Ax) - D<j>(x) ■ n(x + Ax)} < \\Ax - Ax\\M*[4>\. 

Therefore it is sufficient to find an upper bound for 



max mm 

p,r A£ 



(D</>(x) -p) ■ Ax + M*[<f>]\\Ax- Ax\\ + -{D 2 (f>(x)Ax, Ax) - ^(TAx,Ax) - e 2 f{t,x, z,p,T) 



In other words, by taking advantage of the monotonicity of the operator S e with (|4.26[) . we shall look for 
an upper bound for an approximated operator bounding S e above and very close to it as e — > 0. 

4.2.1. Proof of Lemma \4-9\ We introduce 
(4.27) 

A x (p, T, Ax) := (D(f>(x) -p) ■ Ax+ \\ Ax- Ax\\M*[4>} + ^ (D 2 <j)(x)Ax, Ax) -^{TAx, Ax)~e 2 f{t, x, z,p, T), 

where Ax = projjy(x + Ax) — x. We give the following useful decomposition: 
(4.28) 

i(L> 2 0(x)Ax, Ax) - i(rAx,Ax) = i(L> 2 0(x)Ax, Ax) - ^(D 2 (f>(x)Ax, Ax) + ^((D 2 (j)(x) - T)Ax, Ax), 
which will be used repeatedly in this section. We clearly have by ()4.20j) that 

|| Ax - Aa;||M*[0] + i(D 2 </)(x)Ax, Ax) - i(L> 2 0(x)Ax, Ax) 

< llAx- Axil ( M 



||Ax-Ax|| ^Af e K [0] + i||L> 2 0(x)||(3e 1 - Q -d(x))^ < 0. 



From the previous inequality and (|4.28|l we deduce that for all p, T, Ax constrained by (|2.16|) - (|2.17l) . 

A x {p, T, Ax) < {D<j>{x) -p)-Ax + \({D 2 <j>{x) - T)Ax, Ax) - e 2 f(t, x, z,p, T). 

By monotonicity of the operator S e and by using (2lL Lemma 4.1] to estimate the max min, we conclude 
that 

S e [x,t,z,4>] - 4>{x) < maxmin \(D<f>(x) - p) ■ Ax + \((D 2 4>(x) - T)Ax,Ax) - e 2 f{t,x, z,p,T) 
p,v Ax L 2 

< -e 2 f(t,x,z,D^{x),D 2 <p{x)) +o(e 2 ), 
which gives the desired result. 
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4.2.2. Proof of Proposition \4-10\ case (i) We define the function A£ of Ax associated to the particular 
choice p = p^ t and T = T opt by 



(4.29) Af(Ax) = {D<j>(x) • Ax + \\Ax - Ax\\M*[<)>] + ±(D 2 q>(x)Ax, Ax) - ±(T opt Ax,Ax), 



where Ax = projjy(.x + Ax) — x. Thus, the operator S e can be written in the form 
(4.30) 

1 

opt — yj ' < -'* L ' ~r 2 



S e [x, t, z, (/)] — 4>{x) — max min 

p,r Ai; 



A x b (Ax) + (pM -p)-Ax + -((r opt - T)Ax, Ax) - e 2 f{t, x, z,p, T) 



To compute an upper bound of (|4.30[) . we now introduce two preliminary lemmas. 
Lemma 4.11. Assume that M*[<f)] > 0. Then Af defined by (|4.29p is Ax-bounded by 

(4.31) < sup^(Ai) < he 1 -* - d(x)) (3Mf [0] + 4p 2 ?!.(x)||e 1 - Q ) , 

Ax 2 

where Ax is constrained by (|2.17p . 



Proof. This estimate follows exactly the same lines as for Lemmas I4.6H4.8I The sup is clearly positive by 

M 

J opt 



considering Ax — 0. Then, by plugging the expression of p^ pt in Ab(Ax), we have 



At(Ax) = { - £ 2£l J X) (Ax)! + \\Ax- Ax\\)M*{<j>] 

+ l( £l_Q " J^)p 2 ^))n( A i)i + \(D 2 <p{x)Ax, Ax) - i(r opt A£, Ax). 
Since M^[<j>\ > 0, using the estimates (|4.7|) and (|4.22[) . we obtained the desired estimate. □ 

Lemma 4.12. Let f satisfy (11. 4p and (II. 6p (1 1 . 7[) and assume a, j3, 7 satisfy (|2.13l) - (|2.15p . Let (p e )o<e<i 
and (r e )o <£ <i be two sequences bounded respectively in M. N and § N . Then for any x, t and z, we have 



max mm 

\\p\\<e-' 3 ||A£||<ei- 

\\r\\<e-y 



( Pe -p)-Ax + ~<(r e - T)Ax, Ax) -e 2 f (t, x, z,p, V) 



= -e 2 f(t,x,z,p e ,T e ) + o(e 2 ). 



Moreover, the implicit constant in the error term is uniform as x, t, and z range over a compact subset 
ofTi xlxl. 

Proof. It is a direct adaptation of Lemma 4.1] by distinguishing three cases according to the size of 
\\p e - p\\ and A min (r e — Y). □ 

We can now provide an upper bound on (|4.30p . By Lemma |4.11| Ab is upper bounded independently 
of all possible moves Ax. It follows from (|4. 30(1 that 



S £ [x, t, z, (f\ — 4>(x) < sup .A^ (Ax) + max min 

Ax P< T A * 



(p%L - P )-Ax-\ (T opt Ax, Ax) -e 2 f (t, x,z,p, V) 



opt if) ^ 2 

The consistency Lemma 14.121 provides an estimate of the max min and one obtains 

S e [x,t,z,</>] - <P(x) < su P A x b (Ax)~e 2 f{t,x,z,p^ t ,T opt )+o(e 2 ). 

Ai 

By plugging the upper bound in (|4.3ip of A% in the previous inequality, we obtained the desired result. 
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4.2.3. Proof of Proposition [^TTD| case (ii) It is sufficient to show that for any ||p|| < e 13 and ||r|| < e 1 , 
there exists ||Ax|| < £ 1_Q , determining Ax by (|2.22|l . such that 

(4.32) (D<p(x) -p) ■ Ax + M*[(j>}\\Ax - Ax\\ + ^(D 2 <p(x)Ax, Ax) - -{TAx,Ax) 

- e 2 f{t, x, z,p, T) < -e 2 f(t, x, z, Dct>(x),D 2 cj>{x)) + o(e 2 ), 

with an error estimate o(e 2 ) that is independent of p and T and locally uniform in x, t, z. In view of the 
conditions (|2.15[) and (|4.23p . we can pick /i > and 6 > such that 

(4.33) fi + 7 < 1 — a and [i + jr < 1 + a, 

(4.34) <5 < min(2a, p - (1 - a)). 
Now we consider separately the following three cases: 

(1) \\D<j>(x) -p\\ < e» and X min (D 2 ^(x) - T) > -e & , 

(2) \\Dcj>{x) -p\\ < e» and X min (D 2 0(x) -T)< -e s , 

(3) \\D<j>{x)-p\\>e». 

For case[Tl we choose Ax = 0. By a reasoning similar to Case 1 in the proof of |2ll . Lemma 4.1], we 
obtained the inequality given by (|4.32| . 

For cases [5] and [3] in order to use the decomposition (|4.28[) . we now give a preliminary inequality. By 
the inequality (|4.20p in Lemma [4~71 we have 

3, 



^(D 2 <f)(x)Ax, Ax) - ^(D 2 <f>(x)Ax,Ax) 



< ^^xMAx-AxWe 1 -*, 



which yields with the assumption M*[<f)] < ^\\D 2 cj)(x)\\e 1 a that 

(4.35) M*[<f>]\\Ax- Ax\\ + ^{D 2 <f>(x)Ax,Ax) - -{D 2 cf>{x)Ax, Ax) < y||DV(aO||e 1_a ||A3 - Ax\\. 

By combining the geometric estimate (|4.4[) with the assumption d(x) > e 1 ~ a — e p , we get that the 
left-hand side of (|4.35l) is upper bounded by ^■\\D 2 (j)(x)\\e 1 ~ a+p . By using the decomposition (|4.28|) we 
deduce that it is sufficient to show that there exists ||Ax|| < e 1_Q such that 

{Dcj>{x) -p)-Ax + \{{D 2 cp{x) - T)Ax, Ax) + E\\ D ^^ x )\\ £ ^+p _ £ z f ( t , x , z,p, V) 

2 D 

< -s 2 f(t,x,z,D(t)(x),D 2 (f)(x)). 

For case [21 we choose Ax to be an eigenvector for the minimum eigenvalue A = X m \ n (D 2 <f>(x) — T) of 
norm e l ~ a . Notice that since / is monotone in its last input, we have 

f(t, x, z,p, T) > f(t, x, z, D 2 4>{x) - XI). 

Choosing Ax as announced, and changing the sign if necessary to make (D(f>(x) —p) ■ Ax < 0, we deduce 
that 

(D<t>(x) -p)-Ax + \-{(D 2 4>(x) - T)Ax,Ax) + Hp 2 ^)^ 1 "^ - e 2 f(t,x,z,p,T) 
2 b 

< L 2 - 2a X + ^WD 2 ^)^ 1 -^ - e 2 f(t, x, z,p, D 2 4>{x) - XI). 
2 b 

If -1 < A < -e s then e 2 ~ 2Q A < - £ 2-2a+5 and f(t,x,z,p,D 2 (j)(x) - XI) is bounded. Since e 1 "^ < 
£ 2-2q+5 by (j^ggp for such x W e have 

\e 2 - 2a X + ^-WD^WWs 1 -^? - e 2 f(t, x, z,p, D 2 <t>{x) - XI) < -] £ 2 - 2a+s + 0(e 2 ). 
2 b 4 

In this case, we are done by (|4.34[) . since the right-hand side is < e 2 f(t, x, z, D(j)(x), D 2 (f>(x)) when e is 
small enough. 
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To complete case [2] suppose A < — 1. Then using the growth hypothesis (II . T[) and recalling that p is 
near Dcj)(x) we have 

\e 2 - 2a \ - e 2 f{t, x, z,p, D*<t>{x) - XI) < ~e 2 - 2a \X\ + Ce 2 (l + \X\ r ). 

Now notice that |A| < C(l + ||r||) < Ce" 7 . Since 7 (r - 1) < 2a we have e 2 " 2a |A| > e 2 |A| r . Therefore we 
deduce by (|4~34|) that 

-i E 2 - 2 «|A| +Ce 2 \X\ r + ^-WD^We 1 -^ < -\e 2 - 2a < -e 2 f(t,x, z, D<j>{x), D 2 0(x)), 
I o 4 

when e is sufficiently small. Case [5] is now complete. 

Finally, to treat case [3] we take Ax parallel to Dcf)(x) — p with norm e 1_Q , and with the sign chosen 
such that 

(D(j>(x) -p)-Ax = -e^WD^x) - p\\ < -e 1 ""^. 

17 

By observing that — \\D 2 (l>(x)\\e l ~ a+p < e^^D^x) - p\\, this case follows exactly the sames lines as 
[2f|, Lemma 4.1]. 



4.2.4. Proof of Proposition \J. 1 U\ case (Hi) This proof is quite similar to case (ii) Since this estimate will 
not be needed in the rest of the paper, we just indicate that we need to distinguish three cases according 
to the respective sizes of \\D(f>(x) — p\\ and A m ; n (-D 2 4>{x) — T) with respect to e M and — C\ — e a , where /i 
is defined by (fQ3]l . 



4.2.5. Proof of Proposition \J. 1 0\ case (iv) This case corresponds to the heuristic derivation presented at 
Section [2.2.11 when m < 0. Recalling that p^ pt and r opt are defined by (|4.14[) - (|4.15|) . our task is to show 
that for any ||p|| < and ||r|| < £~ 7 , there exists ||Ax|| < e 1_Q , determining Ax by ()2.22p . such that 

(4.36) {D4>{x) -p) ■ Ax + \\Ax - Ax\\M*[<f>] + ^(D 2 <p(x)Ax, Ax) - - (TAx, Ax) 

-s 2 f(t,x,z,p,T) < -(e 1 — -d(x))M*[4>]~£ 2 min f(t, x, z,p, T opt ) + o(e 2 ), 

with an error estimate o(e 2 ) that is independent of p and T and locally uniform in x, t, z. We can notice 
in (|4.36[) that the function M^[(j>] is e,x-boundcd by ||/i||l°c + |Z?<j!>||l~. Moreover, by Lemma |4~T1 we have 



(4.37) 



l -{D 2 ^{x)Ax,Ax) -±(r opt Ax,Ax) < ^||D a ^(x)||(e 1 -« - d(x))e 1 -". 



Thus, it is sufficient to examine, for any ||p|| < e 13 and ||r|j < e 7 , 



(4.38) min 

Ax 



(P4>(x)-p) ■ Ax + \\Ax- Ax\\M*[<t>] + -((r opt - T)Ax, Ax) - e 2 f(t,x, z,p,T) 



We consider separately the following three cases: 

a. \\pM t -p\\ < 3(l - 0) \M*[<t>]\, and A min (L opt - Y) > -e a , 

b. \\ p v t -p\\ < 3(l - |^r) |M*[0]|, and A min (L opt - L) < -e° 
c \Kt-p\\>2(l-0l)\M?[<t>]\. 

For case (a) we choose Ax — ±e 1 ~ a n(x) with the sign chosen such that 

o 



(p-p™)-Ax<0. 
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Since A m in(r opt — T) > — e a we have r opt — T + e a I > and thus T < T opt + e a I . Using the monotonicity 
of / with respect to its last entry, this gives f(t,x,z,p 7 T) > f(t,x,z,p 7 T op t + e a I), Since / is locally 
Lipschitz, we conclude that 

(4.39) f(t,x,z,p,T)> f(t,x,z,p,T opt )+0{e a )> min f(t, x, z,p, T opt ) + 0{e a ). 

P6-B(p^ t ,r) 



The constant in the error term is independent of p and T, since we are assuming in case (a) that \\p— p^\\ < 
3(||/i||l°° + ||Z>(^)||ioo). Moreover we directly compute 

(4.40) (D4,(x) ■ Ax + \\Ax - Ax\\M*[<f>] = ^(e 1 ^ - d(x))M*[<t>]. 
Since e x ~ a — d(x) > e p and M*[<j>] < 0, we have 

(4.41) -(e^ - d(x))M?[<f>] < ~£ P M£[<I>] < - 1 -s^^+p. 
By noticing that £ 2 - 2q -7 < £ i-«-«+p usmg ([j^j]) , we deduce from (|4~4"Tj) that 

(4.42) |i((r opt -r)Ax,Ai)| < \{\\D 2 <f>{x)\\ +e-i)e 2 - 2a < ^s 2 - 2a ^ « {e 1 ^ - d{x))M^]- 

Therefore, by combining (|4.37|) . (|4.40p and (|4.42p , the choice Ax = ±e 1_Q n(x) in the left-hand side of 
(|43o) yields 

{Dcj){x)-p) ■ Ax + ||Ax- Ax\\M*[<f>] + ^(D 2 cf>(x)Ax, Ax) - ~ (TAx, Ax) — e 2 / (t, x,z,p, T) 

< i(£ 1 - Q -d(x))(AffM + ^||^ 2 0(x)||e 1 -«) + fe 2 - 2Q -''-e 2 mm f(t, x, z,p, T opt ) + o(e 2 ) 

<-(e l - a -d{x))M*[<j>]-e 2 mm f(t, x, z,p, T opt ) + o(e 2 ), 

4 peB(pM t ,r) 

as desired. 

For case |(b)[ in view of the condition (|4.23[) , we can pick a > 1 — a such that 

(4.43) p<(T <l_2il_!2. 

Let v x be a unit eigenvector for the minimum eigenvalue A = A m i n (r opt — T). We choose Ax of the form 

(4.44) Ax = ± [(e 1 "" - e a ) n(x) + sgn((n(x), v x ))e a v x ] = ± [ ai n(x) + bv x ] , 

where a\ — (e 1_Q — e a ), b — sgn((n(x), v x ))e a and sgn denotes the sign function with the convention 
that sgn(0) = 1. The sign ± will be chosen later. This move fulfills the following estimate. 

Lemma 4.13. The move Ax defined by (I4.44|) is authorized by the game and satisfies 

(4.45) (D<f>(x) ~p% t ) ■ Ax+ \\Ax - Ax\\M^} < - - ~ d[x) (M°[fl + p 2 0(x)||e 1 -«) - 4e"M*[0], 
independently of the choice on ± in (|4.44p . 

Proof. To authorize this move, it sufhces to check that ||Ax| < e 1_Q . After some calculations and by 
rearranging the terms, we compute 

||Ax|| 2 = e 2 - 2a + 2e 2CT - 2s 1 ~ a+ ' T + e a (s l - a - e CT ) \(n(x), v x )\ 

= e 2 - 2a - 2e 1 - a+a {l - e ff - 1+a )(l - i|(n(x), < e 2 " 2 ". 

For the second part, we distinguish successively the two cases ±. By (|4.44p . we directly compute 

(4.46) Ax • n{x) = ± [(e 1 -" - e a ) + \(n(x),v x )\e a ] = ± [e 1 "" - (l - |(n(x), v A )|) e a ] . 
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If Ax ■ n(x) < 0, this move corresponds to the sign — in (|4.44|) by (|4.46[) and we observe that x £ f2 by 
Lemma 14.21 As a result, by introducing the explicit expressions of p^ t and (Ax)i respectively given by 
(|4~14| and (|4~46|) , we get 



(447) (D</>(x) -pf pt ) • Ax + \\Ax - Ax\\M*[4>] = {D<t>{x) - < t )i(Az)i 

= - (-5(1 - §^WM + \{e l - a - d ^){D^{x))n) (e 1 - - (1 - \(n(x),v*}\)e° 
Since < e l ~ a - (1 - \(n(x),v x )\)e cr < e 1 "", we observe that 



\(^ a ~ ^M)(D 2 Hx)) n (e 1 - - (1 - Kn(x),« A ) 



2 

By plugging this inequality in (|4.47p and rearranging the terms, we obtain 
(D4,(x) - < t ) • Ax+H Ax - Ax\\M x M 



<±\\D^(x)\\(e 2 - 2a -d 2 (x)) 



< (e 1 "" - d(x)) -(!-(!- \(n(x),v x )\)s^ 1+a ) MM + dl^MHe 1 



< \{z l ~ a d{x)){M*M + WD^WWe 1 -*) \e°M*M- 
Otherwise, if Ax ■ n(x) > 0, this move corresponds to the sign + in (|4.44p by (|446p . We have 

\\Ax-e l - a n{x)\\ = \\-e a n[x) + sga({n(x), v x ))e a v x \\ = V2e a Jl - \(n(x),v x )\ < V2e CT . 

By using Lemma T4.21 we deduce from the previous inequality that, for e small enough, the intermediate 
point x — x + Ax is outside Q and 

(448) e 1 -" - d(x) - V2e a - 2C l£ 2 ~ 2a < \\Ax - Ax\\, 

where C\ is a certain constant depending on the principal curvatures of dQ in a neighborhood of x. By 
repeating the computations above, we find 

{D<t>{x) -p&MAa:)! < ^(e 1 ^ - d(x)) {-(1 - (1 - \(n{x),v x )\)e°- 1+a )M*[<t>] + WD^x)^ 1 -"} 

< \^- a - d(x))(-M*[J>] + WD'^We 1 -*). 
Recalling that [</>] < 0, by combining (|4.48|) with the previous estimate, we are led to 
(Dtftx) ■ Ax+ \\Ax- Ax\\M x M 

< ^(e 1 -" - d(x))(-M*[<f>] + H^^lle 1 -") + (e 1 "" - d(x) - sf2e° - 2C 1 e 2 - 2a ^j M*[4>] 

< ^(s 1 - - d(x))(M?[<l>] + ||£>V(^)l|£ 1_a ) ~ e°M*[4>](V2 + 2C 1 E 2 - 2a - a ). 

Putting together the two cases, the proof of the inequality given by (|4.45|) is complete. □ 



Now we turn back to the analysis of case |(b)[ Note that since / is monotone in its last input 

f(t,x,z,p,T) > f(t,x,z,p,T opt - XI). 
The direct evaluation of the second order terms in Ax of (|4.38p gives 

((r opt - T)Ax , Ax) = al({r opt - T)n(x), n(x)) + 2 ai b((T opt - T)v x ,n(x)) + b 2 ((T opt - T)v x ,v x ) 
< Oi(||ropt|| + ||r||) + 2 ai b\(v\ n(x)) + b 2 X. 
With our choice for Ax, we have aib(v x , n(x)) > 0. Hence, since A < in case |(b)| it follows that 
((r opt - T)Ax, Ax) < a 2 (\\T opt \\ + ||r||) + b 2 X < e 2 - 2a {\\D 2 ^x)\\ + e^) + e 2 °X. 
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Choosing Ax as announced, using (|4.37|) and (I4.45[) and changing the sign ± in (|4.44|) if necessary to 
make (p opt — p) ■ Ax < 0, 

{D(j>{x) -p) ■ Ax + \\Ax — Ax\\M*[4>] + i(Z) 2 0(x)Ax, Ax) - ~{TAx,Ax) - e 2 f (t,x,z,p,T) 

< l( e i-«_ d(a;)) ^A/n0] + ^ 2 ^)lk 1 - a ) +^ 2 - 2Q (ll^ 2 ^)ll +e^)-4 £ CT M^] 
(4.49) + ie CT A - e 2 f(t, x, z,p, T opt - XI). 



Since d(x) < e 1 a — e p in case (iv) we deduce from the assumption (|4.43|) that 

(4.50) e 1 - 01 - d{x) > e p e a . 

Since M*[(j>] < -e 1 -"^ and £ 2 - 2a -~t ^ £ i-a- K + P using we conc i u d e by (j430| that 

(4.51) i( £ 1 -«-d( a; ))(MJ ! [0] + 5|p 2 ^)lk 1 - Q ) + ^"(p 2 ^)!! + £ ^)-4 £ ^H 

<\{e l -« - d{x))M*[<])}. 

It remains to control the terms in (I4.49P depending on A. If — 1 < A < ~e a , then e 2a X < —£ 2a + a and 
/(£, x, z,p, r opt — XI) is bounded. So for such A we have 

(4.52) tj£ 2<t A - e 2 f(t, x, z,p, T opt - XI) < -\e 2a+a + 0(e 2 ). 

In this case, the right-hand side is < —e 2 min f(t,x,z,p,T opt ) when e is sufficiently small since 

Pei?(p« t ,r) 

e 2a+ a > £ 2 by ((^431 

To complete case |(b)| suppose A < — 1. Then using the growth hypothesis (jl.7|) and recalling that p is 
near p opt we have 

(4.53) ie 2ff A - e 2 f(t, x, z,p, D 2 J>{x) - XI) < ~£ 2ct |A| + Ce 2 (l + |A| r ). 

Now notice that |A| < (7(1 + ||r||) < Ce~^. Since 7 (r - 1) < 2 - 2a by (14^31 . we have e 2<T |A| > e 2 |A| r . 
Therefore 

-\e 2 °\X\+Ce 2 \X\ r <~-e 2 ° <-e 2 min /(t, as, z, p, T opt ), 
for e small enough. Case|(b)|is now complete. 



Finally in case (c) we take Ax to be parallel to p^ t — p with norm e 1 Q , and with the sign chosen 
such that 

(4.54) « t -p) ■ Ax = -s^WpZU-pW < -3(e 1 -" - d(x))\M* [</>]\ < -Se 1 -"-^. 

Estimating the other terms on the left-hand side of (I4.36[) , some manipulations analogous to those made 
in Lemma \4. Ill led us to 

\{D4>{x) -p ^J • A*+ ||Ax- Ax||M*[0]| < \{e l ~ a - d{x)) (3|M* [0]| + 4||# 2 0(x)||e 1 -«) . 



From (|4. 54(1 . we deduce that 



_l-a||„M „ll , oil n2ic„Ml J-2a 



(D0(z) -p) • Ax + || Ax - Ax||M*[0] < --e 1 ""!^ - p|| + 2||^ 2 0(x)||e 
Estimating the other terms 

(4.55) K(r opt (x)-r)Ax,Ax)| < (C+||r||)||Ax|| 2 < Ce-i +2 - 2a , 
and 

(4.56) e 2 \f{t, x, z,p,T)\ < Ce 2 (l + ||p||« + \\T\\ r ) < C{e 2 + e 2 ||p||« + e 2 ^)- 
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Thus 

(Dcj)(x) -p) ■ Ax + \\Ax - Ax\\M*[(j)} + ^(D 2 cj)(x)Ax 7 Ax) - UvAx, Ax) - e 2 f (t, x, z,p, T) 

< -^ 1_Q bopt - P\\ + Ce 2 |b|| 9 + 0{e 2 - 2a + £ -~i +2 - 2a + e 2 -~> r ). 
Since e l - a \\p^ pt - p\\ > 2e 1 - a - K+p by using ([4"33)l . we obtain that 

(4.57) £ _ 7+ 2-2a + £ 2- 7 r <<; £ X-a|| p M _ p y ( 

noticing that min(— 7+2— 2a, 2— 77") > 1— a— K+p by using (|2.15l) and (|4.24p . Thus, by combining (|4.55[) - 
(|4.57p . we conclude that 

{D(j){x) -p) ■ Ax + || Ax - Ax\\M*[(j)] + ^(D 2 (/)(x)Ax,Ax) - ~(TAx,Ax) - e 2 / (t,x,z,p,T) 



< -- 



1 



^Ib^t-pll + C^lbH 



2V2 

If ||p|| < 2||p^ t ||, then e 2 |b|| 9 < £ i-«-«+P. if || p || > 2||p^ t ||, we infer from the condition on /3 in ([2~T5)l 
that e 1_a ||pQp t — p\\ ~ £ 1_Q ||p|| ^ £ 2 ||p|| 9 . In either case the term e 1 ~ a \\p^p t — p|| dominates and we get 

(p£t -p) • Ax + i((r^ t - r)Ax, Ax) - e 2 f(t,x,z,p,r) < -\e l -<*\\p™ t -p\\ < -^{e 1 -* - d{x))M*[<t>). 

The right-hand side of this inequality is certainly < -(e 1_Q — d(x))M* [tp] — e 2 min /(£, x, z,p, r opt ) 
when e is small. Case (c) is now complete which finishes the proof of Proposition 14.101 

4.3. Application to stability. To prove stability in Section [5l we will need some global variants of 
Propositions 14.51 and 14.101 It is at this point that the uniformity of the constants in (|1.6j) - ()1.7j) in x and 
t, and the growth condition (|1.7[) intervene. We must also take care of the Neumann boundary condition. 
Unlike the Dirichlet problem solved in [2l(, it is no longer appropriate to consider constant functions as 
test functions. For this reason, we are going to consider a C 2 (£!)-function tp such that 

(4.58) = \\h\\ L oo + 1 ondn. 

on 

It is worth noticing that ip has exactly the same properties as the function introduced in Section ^. 3l for the 
game associated to the elliptic PDE with Neumann boundary condition. If we take tp = (||/i||l°° + l)ipi 

dtp 

where ipi <E C 2 (f2) such that — — = 1 on dQ,, it is clear that HV'llc^TT) = ll^ 1 llc 2 (o) (•"■ + ll^lk 00 )- 
The next lemma is the crucial point to obtain stability in both parabolic and elliptic settings. 

Lemma 4.14. If ' %p G C 2 (Q) satisfies (|4.58p . then there exists Sq > such that for all e < Eq and for all 

xefi(e 1 - Q ) ; 

(4.59) -\\ h \\ L «,-\\Di;\\ Loom <M*[i,}<-± and \ < m%[-^\ < \\h\\ L - + ||L>V>IU (n) - 

Proof. We shall demonstrate the bounds on M'£ [ip] in f|4. 59[) : the proof for m%\— tp] is entirely parallel. 

The left-hand side inequality is clear by the Cauchy-Schwarz inequality. Let us consider < s < Eq, 

1 

where Eq = (^\\D 2 tp\\ Lcc ^ + 2^ 1 " . By the geometric relation (|4.4[) . we observe that every move Ax 
associated to the move Ax decided by Mark satisfies 

IIAxll < 2e 1 - a < ' 



2\m\\ L °°(n } + 1 
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By the Cauchy-Schwarz inequality and using that tp € C 2 (Cl), we have 

h(x + Ax) - Dij}(x) ■ n(x + Ax) < \\h\\L°° - D^(x + Ax) ■ n(x + Ax) + (Dip(x + Ax) - Dip(x)) ■ n(x + Ax) 

<-l + \\D 2 i,\\ Loo(n) \\Ax\\<-~. 

Then, by passing to the sup, we get the desired result. □ 

Lemma 4.15. Let <p E C 2 (£l). Assume that p™ pt p^f pt and T opt are the strategies, associated to (j>, 
respectively defined by P~T5|l . P~14]) and (|4~T3]l . Then, for all x E n(e 1 ^ a ), we have 

3, 



max(||p* t (i)|U|j^(x)||) <i(||fe|| L -+3||D^|| c i (n) ) and \\T opt (x)\\ 



< -WD 2 



Proof. The proof being exactly the same for p™ t , it is sufhcient to show the result for p™ t . By the 
triangle inequality and (|4.14p , we have 



\\p% t (x) - D<p(x)\\ < i (l - (\MM\ + Is 1 - (l + \\D 2 <f>(x)\\ 

K^dM^m+e^WD^x)]]). 

Since M^[(j>] is e,a;-bounded by ||/i||l°° + I-D^IIl^TTp we deduce the desired inequality on ||p^t(aj)||. 
Similarly, the estimate on ||r pt(x)|| stems directly from (|4.15| and the triangle inequality. □ 

In preparation for stability, we need to compute the action of S e on ip. According to Lemma |4.14| only 
some cases proposed in Proposition 14.101 must be considered. The next proposition gives the required 
estimates for S e concerning these cases. 

Proposition 4.16. Let f satisfy (TO]) and ([TB] ) -([T7 )l and assume a, f3, 7, p fulfill ([235]) - ([235]) and 
(|4.23l) . Then for any x, t, z and any C 2 (fl) -function (j) defined near x, S e [x, t, z, <j)\ being defined by (|4.1I1 . 
we have 

(4.60) S e [x,t,z,4>]-<t>(x) 

'Ce 2 (l + \z\), ifd{x)>e 1 - a , 



< < 



3s 1 - a M^] + Cs 2 (l + \z\), ifd(x) < e 1 "" and M*[4>] > |||D 2 ^(x)||e 1 - a ) 
Ce 2 (l + \z\), ife 1 -* - e" < d(x) < e 1 ^ and M*[(/>] < |||Z> a ^(a;)||e 1 - 01 , 

Mf>M%[<l>] + Ce 2 (l + \z\), ifd{x) < e Y - a - s p and M*[<f>] < -e 1 - Q " K , 



with a constant C that depends on + but is independent of x, t and z. 

Moreover, if d(x) > e^ a , or if d(x) < e^ a and m%[<j>] > ^(Se 1 -" - d(x))\\D 2 <j)(x)\\, then 

(4.61) -Ce 2 (l + \z\) < S e [x,t,z,<j)]-<l>{x), 

with a constant C that depends on \\D<j)\[ c ifj^ but is independent of x, t and z. 

Proof. The arguments in the different cases are the same as those given in the proof of Proposition 14.101 
but we must pay attention to the uniformity of the constant. For the second part, since / grows linearly 
by CU) and \\(D(f>(x), D 2 (p(x))\\ < \\D<t>\\ cm , we have 

(4.62) \f(t,x,z,D0(x),D 2 0(x))\ <C(l + \z\), 

with a constant C that depends on ||£)0|| c i^ but is independent of x, t and z. The lower bound 

S E [x, t, z, 0] ~ 4>(x) > -e 2 f{x, t, z, D(/)(x),D 2 <j)(x)) > -Ce 2 (l + \z\) 
is a consequence of Proposition 14.51 and (|4.62[) . 
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Similarly, since we know by Lemma 14.151 that max(||p™ t (x)||, HPoptl 2 -)!!) + ||r op t(^)|| is uniformly 
bounded by gll^lli 00 + 3||-D<^||c?i(Q), we get that 

(4.63) max(|/(t,a:,«,p^(i) j r < ^(s))U/(i,i,2,p^ t (i),ro pt (a:))|) <C(l + \z\), 
with a constant C that depends on H-D^H^i^ and ||ft||i«> but is independent of x, t and z. 

We shall prove the estimate for the fourth alternative of (|4.60|) by examining the proof of Proposition 
14.101 case (iv) the proofs for the other alternatives being quite similar. Since / is locally Lipschitz by 

f(t,x,z,p,T ovt )>f{t,x,z,pM t {x),T ovt {x))-C{l + \z\) + 11^11^(0)) , 

where C depends only on 1 1 X) e^> 1 1 ^ ! (T2 ) anc ^ II^IU 00 by the estimates on p^ t and r opt given by Lemma l4.15l 
By using (|4.63l) . we deduce that there exists a constant C depending only on ||Z?0|| c i^ and such 
that 

(4.64) min f(t, x, z,p, r opt ) > -C(l + \z\). 



mm 

PGB(p^ t ,r) 



In case (a) by combining (|4.64[) and the locally Lipschitz character (II. 6p of / on T, the estimate (|4.39D 
gets replaced by 

f(t,x,z,p,T)>-C(l + \z\)(l+e a ), 
whence by (|4.63[) there exists a constant C depending on ||-D<^|| c i^ + such that 

-e 2 f(t,x,z,p,T)<C(l + \z\)e 2 . 

In case |(b)| since the domain satisfies both the uniform interior and exterior ball conditions, we notice 
that the constant C\ corresponding to the curvature of the boundary (see Lemma I4.2[) is x-bounded. 
This implies that the first order estimate (|4.45p is valid independently of x for e sufficiently small. Thus, 
the estimate (|4.51j) is valid uniformly in x. Besides, the estimate (|4.52p gets replaced by 

i £ 2CT A - e 2 f(t,x, z,p, T opt (x) - XI) < -\e 2a+a + Ce 2 (l + \z\)\\p\\ \\T opt (x) - XI\\, 

where C depends on ||-D^||c?i(fj) + II^IU 00 - We obtain an estimate of the desired form by dropping the 
first term and observing that A is bounded. In second half of case |(b)| and in case (c) we used the growth 
estimate (JT77J) ; since z enters linearly on the right-hand side of (|1.7[) . the previous calculation still applies 
but we get an additional term of the form C\z\e 2 in (|4.53[) (|4.56p . □ 

The following corollary provides the key estimate for stability in the parabolic setting. 

Corollary 4.1. Let f satisfy (U]) and (fTTB|) - ([TT7|) and assume a, (3, 7 fulfill ([2~T3|) - (|2~l"5l . Then, for 
any x, t, z and ip 6 C 2 (f2) satisfying (|4.58| . we have 

(4.65) S £ [x,t,z,ip] -if)(x) < C(l + |z|)e 2 and S e [x, t, z, -1/3] - (-ip)(x) > -C(l + \z\)e 2 , 
with a constant C that is independent of x, t, z but depends on H-DV'llc^n) an< ^ ll^-lli 00 - 

Proof. We shall prove the first estimate, the second follows exactly the same lines. By applying Lemma l4.14[ 
we have that M*[ip] < — | for all x <G r2(e 1_a ). We introduce p fulfilling (|4.23[) . By putting together the 
estimates obtained from (|4.61[) and the third alternative in (|4.60p . we get that there exists a constant C 
depending only on ||-D^||ci(Q) and such that 



S £ [x, t, z, ip] — ip(x) < 



Ce 2 (l + \z\), ifd(x)>e 1 - a 
\ePM^] + Ce 2 (l + \z\), if d(x) < e 1 ^ 



Noticing that M*[%l)] is negative, we get the proposed result. □ 
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4.4. The elliptic case. For the game corresponding to the stationary equation, we consider the operator 
Q £ defined for any x £ Q, and any continuous function <f>: — > R, by 

(4.66) Q e [x,z,<f>] =supinf \e~ Xe2 <f>(x + Ax) 

- (p- Ax + i (TAx, Ax) + e 2 f{x,z,p,T) - ||Ax - Ax\\h(x + Ax) 

with the usual conventions that p, T and Ax are constrained by (12.161) and (I2.17P and that Aa; is 
determined by (|2.22[) . We can easily check that the operator Q e is still monotone but its action on shifted 
functions by a constant is described by the following way: for all function <f> E C(£l) and cel. 

(4.67) Q e [x, z, c+(j>] = e- Xe \+ Q £ [x, z, 4>] . 

The dynamic programming inequalities (|2.6ip - (|2.62p can be concisely written as 

u s (x) < Q e [x, u e (x), u e \ and v £ (x) > Q s [x, v £ (x), v 6 ]. 

In the elliptic setting, we can formally derive the PDE by following the same lines as for the parabolic 
framework. We keep the optimal strategies p^ t , p^t and r opt for Helen, defined by (|4.13[) . (|4.14[) and 
(I4.15P in an orthonormal basis 23 = {e\ = n(x), e2, • • • , ejv). The next proposition is the elliptic analogue 
of Propositions 14.51 and 14.101 It establishes the consistency estimates for Q £ defined by (|4.66[) . 

Proposition 4.17. Let f satisfy (|1.4p and (|1.9p - (jl.lOp and assume a, j3, 7 and p fulfill (|2.13p - (|2.15p and 
(I4.23p . Let p™ pt , p^ pt and T opt be respectively defined in the orthonormal basis 25 = (e\ = n(x), e2, • • • , eAr) 
by (|4.13p - (|4.15p . For any x, z and any smooth function <fi defined near x, we distinguish two cases for 
the lower bound estimate: 

i. Big bonus: if d(x) > e x ~ a or m x e [(f>] > \(is 1 ~ a — d{x))\\D 2 cj)(x)\\, then 

(4.68) - e 2 {f{x, z, D(f>(x),D 2 (f)(x)) + \<f>(x)) < Q £ [x, z, 0] - <j>{x). 

ii. Penalty or small bonus: if d(x) < e 1 ~ a and m*[(f)] < i(3e 1_Q — d{x))[\D 2 (f){x)[\, then 

^(e 1 -" - d(x)) (sm x M - 4|| J D 2 ( Xx)|| £ 1 -") - e 2 (f(x, z,p™ t (x), T opt (x)) + \cf>(x)) < Q e [x, z, <f>] - 

where s = — 1 if [(f)] > and s = 3 if [(f)] < 0. 

For the upper bound estimate, we distinguish four cases: 

i. Big bonus: if d(x) < e 1 ^ and M*[<f>] > § \\D 2 <j>(x)\\£ 1 - a , then 

Q e [x, z, 4>] - </>(x) < 3(e 1 - Q - d{x))M*[<t>] - e 2 (f(x, z,p% t (x), T opt (x)) + \<j>(x)) + o(e 2 ). 

ii. Far from the boundary with a small bonus: if e 1 ~ a —e p < d(x) < e 1_a and [(f)] < ^\\D 2 (f>(x)\[£ 1 ~ a , 
or if d(x) > e x ~ a , then 

(4.69) Q e [x, z, (f)] - (f>(x) < -e 2 (f(x, z, D<j>(x), D 2 (f){x)) + X(f>(x)) + o(e 2 ). 

Hi. Close to the boundary with a small bonus /penalty: if d(x) < e 1 ^" — e p and — e 1 ~ a ~ K < M^[<p] < 
IIID 2 ^^)^ 1 -", then 

Q e [x, z, 0] - <f>{x) < -e 2 (f(x, z, Dcf>(x), D 2 (f>{x) + C X I) + \<j>(x)) + o(e 2 ), 

with & = f i^o^ii 

iv. Close to the boundary with a big bonus: if d(x) < e x ~ a — £ p and M^[(f)] < — e 1 ^ a ^ K , then 

(4.70) Q e [x,z,4>]-<j>(x) < i(e 1 - a -d(x))M £ a; [0]-£ 2 ( min f(x, z,p, T opt (x)) + \cj>(x) J + o(e 2 ), 

4 - \peB(pM(x),r) 



with r defined by r = 3 (l - \M*[<j>}\. 
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Moreover the implicit constants in the error term are uniform as x and z range over a compact subset of 
UxR. 

Proof. The arguments are entirely parallel to the proofs of Propositions 14.51 and 14.101 □ 

For stability we will need a variant of the preceding lemma. This is where we use the hypothesis (jl.8p 
on the z-dependence of /. 

Lemma 4.18. Let f satisfy (|1.4| and (jl.9| - (|1.10ll and assume as always that a, (3, 7 satisfy (|2. 13|) — 
(|2.15p . Letip £ C 2 (fl) satisfy (|2.55p . Fix M andm two positive constants such that m+2\\ip\\ Loa ^ < M. 

Then, there exists = C»(||-D'f/'llc' 1 (n)> 11^11^°°) sucn that for any \z\ < M and any x £ CI, we have 
Q E [x,z,m + ip] - (m + ip(x)) < e 2 (1 + (A — rf)\z\ + C*) - Ae 2 (m + ip(x)) , 

and 

Q e [x, z, -m - V] - (-m - ip(x)) > -e 2 (1 + (A - rj)\z\ + C*) - Ae 2 (-m - ip(xj) , 
for all sufficiently small e (the smallness condition on e depends on M, but not on x). 

Moreover, if <j> £ C' 2 (fl), then there exists C = C(M, ||-D<^||ci(f})j ll^lli- 00 ) such that for any \z\ < M 
and any x £ ft such that d(x) < e 1 ~ a - e p and M£[<f>] < -e 1 ~ a ~ K , 

(4.71) Q s [x,z,cf>] - 4>{x) < ~ (e 1 -" - d(x)) M*{4>] + Ce 2 - \e 2 4>(x), 
for all sufficiently small e (the smallness condition on e depends on M, but not on x). 

Proof. We shall prove the first inequality, the proof of the second being entirely parallel. The assumption 
\z\ < M ensures that the constants in (|1.9[) and (ll.lOj) arc uniform. Then the implicit constants in the 
error terms of (|4.69p and (|4.70p are cc,z-uniform for e small enough, and the smallness condition depends 
only on M. Since m + 2||?/>|| Loo( jj) < M we can use the dynamic programming inequalities (|2.61| - (12.62[) . 
First of all, by the action of Q E on constant functions provided by (|4.67l) , we have 

\ 2 

Q e [x, z,m + ip] — (m + ip(x)) = (e~ 6 - l)m + Q e [x, z, tp] — ip(x), 

and noticing that e~ Xe m = (1 — Ae 2 )™ + 0(e 4 m), it is sufficient to get the estimate corresponding to 
rn = 0. By Lemma [4.14[ we observe that every x £ fl(e 1 ~ a ) satisfies M^[tjj] < —i. We now need to 
distinguish two cases according to the distance to the boundary by introducing p fulfilling (|4.23p . If x £ fl 
such that d(x) > e 1- " — e p , since |j (Di/)(x), D 2 ip(x))\\ < K\ — ||-DV'llc 1 (?2)' we deduce by assumption (|1.8p 
on / that there exists such that for all x we have 

\f(x,z,D^(x),D 2 ^(x))\ < (\- v )\z\+C* Ki , 

which gives by (|4.69p that for all x £ fl such that d(x) > 

(4.72) Q e [x, z, V] - ip(x) < e 2 ((A - n)\z\ + C* Kl ) - Xe 2 ^(x) + o(e 2 ). 

If x £ fi such that d(x) < e x ~ a — e p , combining the triangle inequality with the inequalities given by 
Lemma 03 gives that, for all p £ B (p^ pt (x), r) with r = 3 (l - |^ \M*[ip}\, 

HO, r opt (x))|| < \\p™ t (x)\\ L ~ +r+ \\T opt (x)\\ L ~ <K 2 = 7 -\\h\\ L ~+6\\DiP\\ cm , 

since M^[ip] is £,£-bounded by + The assumption (| 1 .8|) on / yields that there exists 

such that, 



(4.73) 



min f(x,z,p,T opt (x)) 

peB(p™ t (x),r) 



<(\- V )\z\+C* K2 , 
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By using this inequality in (|4.T0[) and recalling that M*[ij)] < —A, we conclude that, for all x € O such 
that d(x) < e x - a - e p , 

(4.74) Q s [x, z, V] - rl>{x) < e 2 ((A - n)\z\ + C£ 2 ) - \e 2 ^(x) + o(e 2 ). 

By comparing (|4.72p and (|4.74p we get the desired result by taking C* = max(C^- i , C* K2 ) . 

To prove the third inequality, it is sufficient to replace the assumption (jl.8[) by (jl.lOp in the previous 
estimates. For instance, instead of (|4.73[) . there exists a constant C depending only on M, ||/i||l°°, and 

f(x,z,p,T opt (x)) 



such that 



mm 

P GB(p"(x),r)" 



< C. The rest of the proof remains unchanged. □ 



5. Stability 

In the time-dependent setting, we showed in Section 13.21 that if v £ and u £ remain bounded as e — > 
then v is a supersolution and it is a subsolution. The argument was local, using mainly the consistency of 
the game as a numerical scheme. It remains to prove that v £ and u £ are indeed bounded; this is achieved 
in Section [5Jl 

For the stationary setting, we must do more. Even the existence of U £ (x, z) remains to be proved. We 
also need to show that the associated functions u £ and v £ are bounded, away from M, so that we can 
apply the dynamic programming inequalities at each x € fi. These goals will be achieved in Section 15.21 
provided the parameters M and m satisfy (i) m = M — 1 — 2||-i/>||z,°° and (ii) M is sufficiently large. We 
also show in Section l5~2l that if / is a nondecreasing function on z then U £ is strictly decreasing on z. As 
a consequence, this result implies that v < u, allowing us to conclude that v = u is the unique viscosity 
solution if the boundary value problem has a comparison principle. 

5.1. The parabolic case. To obtain stability, we are going to consider one more time a C^(f2)-function 
ip such that 2^ = + 1 in order to take care of the Neumann boundary condition. 

Proposition 5.1. Assume the hypotheses of Propositions \4-5\ and \4-10\ hold, and suppose furthermore 
that the final-time data are uniformly bounded: 

\g(x)\ < B for all x e H. 
Then there exists a constant s = s(|| ^ He 2 ($?))> independent of E, such that 

u £ (x,t) < (B+ |Mli°a(n))s T ~' + ip{x) for all ie!!, 

and 

v £ {x,t) > —{B + ||V'llL=o ( n))s T " t - tp(x) for all icH, 

for every t < T . 

Proof. We shall demonstrate the lower bound on v e \ the proof of the upper bound on u e is entirely 
parallel. The argument proceeds backward in time tu = T — ke 2 . At k = 0, we have a uniform bound 
v e (x,T) = g(x) > — B by hypothesis, and we may assume without loss of generality that B > 1. Since tj) 
is bounded on Q, we can suppose that 

v E (x,T) =g(x) > -B -il?(x), 

where Bq = B +\\ijj\\ Lao Now suppose that for fixed k > we already know a bound v £ (-, tk) > —B^—ip. 
By the dynamic programming inequality (|2.25p . we have 

v £ (x,t k ~e 2 ) > S £ [x,t,v £ (x,t k -e 2 ),v £ (.,t k )] . 

Since S e is monotone in its last argument, we have 

v £ (x, t k - e 2 ) > S e [x, t, v £ (x, t k - e 2 ), -B k - tp) ■ 
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By applying successively (|4.2[) and Corollary 14. 11 we deduce that 

S E [x,t,v £ (x,t k -e 2 ),-B k -if\ = -B k + S £ [x, t, v £ (x, t k ~ e 2 ), -ip] 

> -B k - i/;(x) - C(l + \v £ (x, t k - e 2 )\)e 2 , 

where C depends only on ||Z3^|| c i^. If v £ {x,t k — e 2 ) > 0, then it is over (recall we are looking for a 
lower bound —B k+ i < —1). Otherwise, we have 

(1 - Ce 2 )v £ (x, t k - e 2 ) > -B k - Ce 2 - ^{x). 

By dividing by 1 — Ce 2 , we get 

B k + C{l + U\\ Lxm )e 2 

Then, by setting B k+ i = — — - — , we obtain 

1 — Ce 2 

v £ (x,t k -e 2 ) > -B k+1 -i/j(x). 
l + C(l + U\\ Lao(Tl) )s 2 

As it is clear that B k+ i < B k —5 — — , we deduce that v e {x, T — fee ) > B k — ip(x) for all 

1 — Ce 2 

k with 

(l + C(l + M Loam )e 2 \ k 
£>fc — -Do — 



1-C*e 2 



Since k — (T — t)/e 2 and recalling that Bq = B + \\ip\\i,oo(fi), we have shown that 

v e {x,t) > -{B+UW^^s^ -ij{x) 

with 

Se ~{ l-Ce 2 ) • 

Since s e has a finite limit as e — > we obtain a bound on v £ of the desired form. □ 

Remark 5.2. By following the construction of the elliptic game we can take ip = (||/i||z,=° + l)^i where 
ipi is defined by (|2.66p . In that case, \\Dip\\ c i^ = + ||/i||z,°°)- This expression can be 

compared for a C 2,a -domain to the estimate given by Remark \2.8\ vrovided by the Schauder theory for 
which WDiplWc 1 ^) pl a V s the role of the constant Cn depending only on the domain. 

5.2. The elliptic case. We shall assume throughout this section that the parameters M and m control- 
ling the termination of the game are related by m = M — 1 — 2||'0|| Loo ^; in addition, we need to assume 
M is sufficiently large. Our plan is to show, using a fixed point argument, the existence of a function 
U £ (x, z) (defined for all x £ CI and \z\ < M) satisfying (|2.58[) and also 

(5.1) -z-x{x) < U £ (x,z) < -z + x(x). 

This implies that U e {x, z) < when z > x( x )i an d U £ (x,z) > when z < ~x(x). Recalling the 
definitions of u £ and v £ , it follows from (12.59p - (12.60p that 

(5.2) \v £ (x)\< X (x), \u £ {x)\< X (x), 

for all x € 0. It is convenient to work with V £ (x, z) = U £ (x, z) + z rather than U £ , since this turns t|5 . 1 [) 
into 

\V £ (x,z)\ < x (x), 



52 



JEAN-PAUL DANIEL 



(5.3) R e [x, z,fa\ = supinf < 



whose right-hand side is not constant. The dynamic programming principle (|2.58[) for II s is equivalent 
(after a bit of manipulation) to the statement that V E is a fixed point of the mapping fa[-, ■) t-¥ R e [-, •, fa 
where the operator R e is defined for any L°°-function <f> defined on x (—M,M) by 

'e- Xe2 c/)(x',z')-S, if \z'\<M, 
- X (x), ifV>M, 
[X(x), iSz'K-M. 

where x' = x + Ax and z' — e Ae (z + 5), with 8 defined as in (|2.57j) . Here p, T and Ax are constrained 
as usual by f2~T6|) - ([2~T7|) . 

We shall identify V s as the unique fixed point of the mapping fa-,-) l— ► ■> ' ■> 4>] m Fx defined by 

(5.4) F x ^{<peL°°(Tlx (-M, M)) :V(x,z) eSJx (-M, M), \fax, z)\ < X (x)} ■ 

Lemma 5.3. Let f satisfy (]1.4p and (|1.9|) - (|1.10p and assume as always that a, f3, 7 fulfill (|2.13p - (|2.15p 
and that fl is a C 2 -domain satisfying both the uniform interior and exterior ball conditions. Then, there 
exists Mq > such that for all two positive constants m and M > Mq satisfying m + 2\\ip\\ LOO (fh = M — 1, 
for any \z\ < M and any x € Q, we have 

Qe[x, z, x] < x(x) and Q e [x, z, ~x] > ~ X (x). 
Proof. We are going to establish the upper estimate for X - By Lemma T4. 181 we deduce that 

Q e [x, z, X ] - X(x) < e 2 (l + (A - rj)\z\ + C.) - Ae 2 (m + ||^|| L „ (n) + i>{x)). 
Since m + 2|j-0|| Loo ^ = M — 1 and \z\ < M, we compute 

Q e [x, z, X ] ~ X(x) < e 2 (l + (A - n)M + C*) - Ae 2 (m - 1 - M L ^ m + iftxj) . 
By rearranging the terms, we obtain that 

Q e [x, z,x] - Xix) < e 2 (l + A(l + 2||V|| i . (n )) + - t]M) . 
We can choose M large enough such that the right-hand side is negative. It suffices to take 

M > M := i (l + A(l + 2||^|| £ao(n) ) + C.) . 
The case for Q s [x, z, — X ] > — x(^) is analogous. □ 



Proposition 5.4. Assume the hypotheses of Lemma \ 5.S\ hold. Suppose further that m — M — 1 — 
211^11 Lao(jiy Then for all sufficiently small e, the map 4>{-,-) <-> R e [-,-,ff>] is a contraction in the L°°- 
norm, which preserves F x . In particular, it has a unique fixed point, which has L°° -norm at most 
m + 2||V>|| i0 c(ny 

Proof. By the arguments already used in (2lL Proposition 5.2], the map is a contraction for any e (this 
part of the proof works for any M). More precisely, if fa, i = 1,2 are two L°°-functions defined on 
Hx (-M, M) to R, then \\R £ [,-,fa]-R e [,-,fa]\\ L - < e'^Wfa - fa\\ L oo. 

Now we prove that if M is large enough and m + 2\\ip\\ L00 ^ = M — 1, the map preserves the ball F x 
defined by (|5 .4[) . Since R £ [x, z, fa is monotone in its last argument, it suffices to show that 

(5.5) R £ [x,z, X ] < X (x) and R e [x, z, ~ X ] > - X (x). 
For the first inequality of (|5.5[) , let p and T be fixed, and consider 

( e- Xs2 X (x') - 5, if|z'|<Af, 

(5.6) in| J - X (x), if z' > M, 

X \ x(x), itz' < —M. 
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If a minimizing sequence uses the second or third alternative then the inf is less than xi x )- I n t ne 
remaining case, when all minimizing sequences use the first alternative, we apply Lemma l5.3l to see that 
(|5.6| is bounded above by x( x )- It follows that for all x € f2, R £ [x,z,x] < x( x )> as asserted. 

For the second inequality of (|5.5|) . the argument is strictly parallel by considering the function — %. 
We have shown that the map (f>(-,-) i-> i? e [-, - ,4>] preserves the ball F x . Since it is also a contraction, the 
map has a unique fixed point. □ 

This result justify the discussion of the stationary case given in Section [5J by showing that the value 
functions u £ and v e are well-defined, and bounded independently of e, and they satisfy the dynamic 
programming inequalities: 

Proposition 5.5. Suppose f satisfies (|1.4|) and (|1.8|) - ||1.9|l . the C 2 -domain fl fulfills both the uniform 
interior and exterior ball conditions, and the boundary condition h is continuous, uniformly bounded. 
Assume the parameters of the game a, /3,7 fulfill (|2.13p - (|2.15[) . ip G C 2 (fi) satisfy (|2.55[) . M /arge enough, 
m — M — 1 — 211-011^00^); X € (0) is defined by (|2.56| . Le£ y £ &e t/ie solution of (|5.3| obtained by 
ProposiUon \5.4\ and let U £ (x,z) — V £ (x, z) — z. Then the associated functions u e , v £ defined by (|2.59p - 
(|2.60p satisfy \u e \ < x an d \v £ \ < X f or a ^ sufficiently small e, and they satisfy the dynamic programming 
inequalities (I2.6ip and (|2.62p at all points x £ O. 



Proof. The bounds on u £ and w e were demonstrated in (|5.2p . The bounds assure that the dynamic 
programming inequalities hold for all i G H, as a consequence of Proposition 12.61 □ 



We close this section with the stationary analogue of Lemma [ 
Lemma 5.6. Under the hypotheses of Proposition HO} suppose in addition that 

f(x, Z\,p, r) > f(x, Zq,P, r) whenever Z\ > zq. 

Then U e satisfies 

U £ (x, z\) < U £ (x, zq) — (zi — zq) whenever z\ > zq. 
In particular, U e is strictly decreasing in z and v £ = u e . 



Proof. The Dirichlet case is provided in 2JJ, Lemma 5.4]. For our game, it suffices to add — ||Ax — 
Ax\\h(x + Ax) in the expression of 5q and 5i defined in the proof of [2l[ Lemma 5.4]. Then the arguments 
can be repeated on the operator R e defined by (|5.3[) . noticing that the function x is independent of z. □ 

6. Some natural generalizations 

In the precedent sections, we solved the Neumann boundary problem in both parabolic and elliptic 
settings. In the present section, we are going to explain without full proof how the previous work can be 
used to solve on the one hand the mixed Dirichlet-Neumann boundary conditions in the elliptic framework 
and on the other hand the oblique problem in the parabolic setting. For the definitions of the viscosity 
solutions on these frameworks which are the natural extensions of those presented in Section I3.1[ the 
interested reader is referred to 12 1 or 



6.1. Elliptic PDE with mixed Dirichlet-Neumann boundary conditions. We extend the games 
of Section 12.31 devoted to the single Neumann problem to the mixed Dirichlet-Neumann boundary- value 
problem 

f(x,u,Du,D 2 u) + Xu~0, in fl, 

(6.1) \u = g, on To, 

du 

— = h, on Tjv, 
on 
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where il C M. N is a domain, Tc U T^v = dfl is a partition of <9f2 with nonempty and closed and 
Tat is assumed to be C 2 . Then, f2 is assumed to satisfy the uniform exterior ball condition and, in 
a neighborhood of Tat, the uniform interior ball condition explained in Definition 11.21 We will need a 
C 2 (f2)-function ip such that 

(6.2) ty = \\h\\ L oo + l onT N . 
From m and ip, we construct a function \ defined by 

(6.3) x( x ) = m + IIV'IU^ + ip(x)- 

As in Section |2"UI we introduce U £ (x, z), the optimal worst-case present value of Helen's wealth if the 
initial stock is x and her initial wealth is —z. The definition of U £ (x, z) for x 6 U Tjy involves here a 
game similar to that of Section l2~3l The rules are as follows: 

(1) Initially, at time to — 0, the stock price is xq — x and Helen's debt is zq = z. 

(2) Suppose, at time tj — je 2 , the stock price is Xj and Helen's debt is Zj with \zj\ < M. Then 
Helen chooses pj g and Tj £ § N , restricted in magnitude by (|2.16|) . Knowing these choices, 
Mark determines the next stock price Xj+i = Xj + Ax so as to degrade Helen's outcome. Mark 
chooses an intermediate point Xj+i = Xj + Axj e R N such that ||A%|| < s 1 ^ 01 . This position 
Xj+i determines the new position Xj+i = Xj + Axj by 

xj+i = proj ?T (i i+ i) G il. 

Helen experiences a loss at time tj of 

(6.4) Sj = pj ■ Axj + - (TjAxjtAxj) + e 2 f(x J ,z j , Pj ,T J ) - \\A Xj - Axj\\h{xj + Axj). 

As a consequence, her time tj+x = tj + e 2 debt becomes Zj+i = e Xs (zj + Sj). 

(3) If Zj + i > M, the the game terminates, and Helen pays a "termination-by-large-debt penalty" 
worth e Xe (x( x j) ~~ Sj) a t time tj+i. Similarly, if Zj+i < —M, the the game terminates, and 
Helen receives a "termination- by-large- wealth bonus" worth e Xe (x( x j) + Sj) at time tj+i- If the 
game ends this way, we call tj + \ the "ending index" tx- 

(4) If 1 2^+1 1 < M and Xj+x 6 Td, then the game terminates, and Helen gets an "exit payoff" worth 
g(xj + i) at time fj+i. If the game ends this way, we call tj+\ the "exit index" ie. 

(5) If the game has not terminated then Helen and Mark repeat this procedure at time tj+i = tj +e 2 . 
If the game never stops, the "ending index" tpc is +oo. 

All the possibilities, apart the end by exit, had already been investigated at Section [2~3l If the game 
ends by exit at time t e , then the present value of her income is 

U E (x ,z ) = -z -S Q - e-^S! e-W-^Se-i + g(x E ) 

= e- XEe \g{x E )- z E ). 

Since the game is stationary, the associated dynamic programming principle is that for \z\ < M, 

~ Xe2 U E (x',z'), iix' £ flUT N and \z'\ < M, 

. -^{g{x')-z'l ifz'el^andlz'l < M, 

(6.5) (7 £ (i,z)=supmm 

pX Ax \- z -x(x), iiz'>M, 

-z + x(x), iiz'<-M, 

where x' — projn(&+ Ax) and z' = e Xe (z + S), with S defined by (I6.4[) . Here p, T and Ax are constrained 
as usual by f2~TC)) - (l^T7|) . 
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The definitions (|2.59[) (|2.60[) of u £ and v £ on U IV are conserved. The corresponding semi-relaxed 
limits are defined for any x £ Q by 

u(x) = lim sup it e (y) and v(x) = liminf v £ (y), 

e^O <=->0 

with the convention that y approaches x from fl U Tjv (since u £ and v £ are only defined on fi U Tjv). 
Proposition 12.61 still holds without any modification for mixed- type Dirichlet-Neumann boundary condi- 
tions. Moreover, the definition of viscosity subsolutions and supersolutions is clear by relaxing the PDE 
condition on T d with the Dirichlet condition in the same way that has been done in pH Section 3] . 

Following the same steps as our proof for the Neumann problem (the main modification consists in the 
proof of convergence on Tjj but has already been done in |2l|). the following theorem is now immediate. 

Theorem 6.1. Consider the stationary boundary value problem (|6.1| where f satisfies (|1.4j) and (|1.8j) - 

(jl.lOp . g and h are continuous, uniformly bounded and Q is a C 2 -domain satisfying the uniform exterior 
ball condition and the uniform interior ball condition in a neighborhood o/Tjv- Assume the parameters 
of the game a, (3, 7 fulfill (pT3]) - ([2T5|l . ip £ C^(U) satisfies ((OJ), \ € C 2 (U) is defined by (Ol), M is 
sufficiently large, and m = M — 1 — 2||T/>|| i=0 ^. Then u £ and v £ are well-defined when e is sufficiently 
small, and they satisfy \u £ \ < \ an d \ v£ \ ^ X- Their relaxed semi-limits u and v are respectively a 
viscosity subsolution and a viscosity supersolution of (|6.ip . If in addition we have v < u and the PDE 
has a comparison principle, then it follows that u £ and v £ converge locally uniformly in f2 to the unique 
viscosity solution of (|6.1|l . 

6.2. Parabolic PDE with an oblique boundary condition. The target of this section is to construct 
a game which could interpret the PDE with an oblique condition h and final-time data g given by 

!d t u - f(t, x, u, Du, D 2 u) = 0, for x £ D, and t < T, 
-^(x,t) =h(x), for x £ dfl and t < T, 

I = g(x), for x £ Q, 

where defines a smooth vector field, say C 2 , on <9f2 pointing outward such that 
(6.7) (s(x), n(x)) >9>0 for all x £ dtl. 

As usual, the domain Q is supposed to be at least of boundary C 2 and to satisfy both the uniform and 
the exterior ball conditions. 



First of all, following P.L. Lions [22], Section 5], P.L. Lions and A.S. Sznitman [2J|, we introduce some 



smooth functions dij(x) = a,ji(x), say C 2 (R N ), such that 
(6.8) 36 > 0,Vx £ R N , (aij(x)) > 9I N , 

N 

\/x £ dft, ajj (x)c,j (x) = ni(x) for 1 < i < N. 
j'=i 

Clearly if we had <; = n, we would just take aij{x) — Sij. Next, the matrices induce a metric d q on R w 
defined by 

1/2 



(6.9) d,{x,y)=mi 



a« (£(«))& (*)€,•(*) 

l<i,j<N 



dt:Z£C x ([o,i},^ N ),m = y,Z(i) = x 



Then it is well known that for ||a; — y\\ small, there exists a unique minimizer in (|6.9p . The interested 



reader is referred to (22| for additional properties about d^. For this specific metric, we can now define 
for any x lying on a small (5-neighborhood of the boundary a unique projection according the vector field 
7 along the boundary by 

(6.10) x 1 = projUV) £ dQ, 
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which corresponds to the unique minimum of d q (x,y) for y lying on the boundary. Finally, B<;(x,r) 
denotes the ball of center x and radius r induced by the metric d s . 

We can now explain the rules of the game corresponding to the oblique problem (|6.6[) . Let the 
parameters a, (3, 7 satisfy (I2.13p - (|2.15p . When the game begins, the position can have any value Xq G CI; 



Helen's initial score is j/o = 0. The rules are as follows: if at time tj 
stock price is Xj, then 



to + J£ 2 Helen's debt is Zj and the 



(1) Helen chooses a vector pj G M. N and a matrix Tj G S N , restricted in magnitude by (|2.16|) 

(2) Taking Helen's choice into account, Mark chooses the stock price Xj+i = a 



Axj so as to 



degrade Helen's outcome. Mark is going to choose an intermediate point Xj+i 
such that 

(6.11) i J+ i€B ( fe,^°), 

which determines the new position + Axj G Cl by the rule 

Xj+i = projl-(x i+ i), 

where proj^ is the projection defined by (|6.10[) . 
(3) Helen's debt is changed to 



Axj G 



pJV 



Zj+l 



(4) The clock steps forward to tj+i — tj + e 2 and the process repeats, stopping when tx = T. At 
the final time Helen receives g{xn) from the option. 

Rather than repeating the arguments already used, we are going to explain the modifications to carry 
out the analysis. First of all, by the boundedness of the and (|6.8|) . the distance d^ defined by (j6.9|) 
is equivalent to the euclidean distance. Since f2 satisfies the uniform exterior ball condition, there exists, 
for a certain r ? > 0, a tubular neighborhood {x G M JV \f2, d(x) < r q } of the boundary on which proj^ is 
well-defined. This guarantees the well-posedness of this game for all e > small enough. Then, if d s or 
the euclidean distance is used to compute D(j> and D 2 <fi for a smooth function <j>, we will get the same 
results. Therefore, we can introduce the oblique analogues £ [<fi] and e [4>] of (|3.ip - (|3.2p by 



(6.12) 
(6.13) 



[</>] := inf {h{x + Ax) - D<f>(x) ■ s(x + Ax)} , 

Ax 

sup {h(x + Ax) - D(f>(x) ■ s(x + Ax)} , 

c+Ax01 
Ax 



Thus, the 



where Ax is constrained by (|6.11[) and Ax is determined by Ax = proj^(x + Ax) — 
particular choices p™ pt , p^p\ and will be now respectively defined in the orthonormal basis 23 ? 
(ei = <;(x 7 ),e 2 , • • • ,e N ) by 



PoptW 

M, / s 
PoptW 



Dcj>(x) 
Dcj)(x) 



; 1- 



M x ) 
d q {x) 



and 



where 



Tl pt (x)=D 2 cj ) (x) + 



-1 



4 



d 2 (x) 



1 - 



d 2 (x) 



d 2 (x) 



{D 2 (t>{x)), 



{D 2 ct>{x)h 



(D 2 4>(x)) 



11 



; [<f>] and M*[4>] are defined by (|6.12[) (|6.13[) . and En denotes the unit-matrix (1,1) in the 
basis 23 s . The definitions of u £ , v £ and their relaxed semi- limits u and v, given by (|2.23p - (|2.24|l and (|2.27|) . 
are conserved. The only change on the dynamic programming inequalities (|2.25l) - (|2.26j) concerning u e 
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and v £ is to replace || Ax — Ax\\ by d^(x + Ax, x + Ax), and to constrain Ai by (|6.1ip . For stability, we 
need to consider a C^(f2)-function tp such that 

^(x) = \\h\\ L ~ + l ondil. 

It is still allowed by the uniform interior ball condition applied to the C 2 -domain Q. By using exactly 
the same ingredients already used for the Neumann problem and adapting the geometric estimates given 
by Section 14.1.11 in the oblique framework, we obtain the following theorem. 

Theorem 6.2. Consider the final-value problem (|6.6p where f satisfies (jl.4p - (|1.7p , g and h are con- 
tinuous, uniformly bounded, fl is a C 2 -domain satisfying both the uniform interior and exterior ball 
conditions, and <; is a continuous vector field on dQ and satisfy (|6.7p . Assume the parameters a, f3, 7 
fulfill (|2.13p - (|2.15p . Then u and v are uniformly bounded on x [i*,T] for any t* < T, and they are 
respectively a viscosity subsolution and a viscosity supersolution of (16. 6p . // the PDE has a comparison 
principle (for uniformly bounded solutions), then it follows that u e and v e converge locally uniformly to 
the unique viscosity solution of (|6.6p . 
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